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Abstract. Two-term Weyl-type asymptotic law for the eigenvalues of one-dimensional 
quasi-relativistic Hamiltonian {—h'^c^d'^/dx'^ + m^c^Y^"^ + Kvou(a;) (the Klein-Gordon 
square-root operator with electrostatic potential) with the infinite square well potential 
Kvcii(a^) is given: the n-th eigenvalue is equal to (n7r/2 — 7r/8)/ic/a -I- 0(l/n). where 2a 
is the width of the potential well. Simplicity of eigenvalues is proved. Some I? and 
L°° properties of eigenfunctions are also studied. Eigenvalues represent energies of a 
'massive particle in the box' quasi-relativistic model. 

O 
(N 

1. Introduction and statement of the results 

^ In quantum physics, the state of a particle, or a system of particles, is represented by 

the wave function \I'(t,x), which satisfies the Schrodinger equation 

th^it,x) = H^it,x). (1) 

(—1 

Q_i Here H is the Hamiltonian, a (typically unbounded) self-adjoint operator, acting in the 

spatial variable x, and h is the reduced Planck constant. When the solution of ([T]) has 
a time-harmonic form \E'(t,x) = e~*'^*/'*$(x), it is called a stationary state with energy 



E. Solutions of this kind are in one-to-one correspondence with eigenfunctions $(x) 

I I and eigenvalues E of the Hamiltonian H. The function $(a;) is often called an energy 

^ eigenstate, and the corresponding E is said to be an energy level of the particle. The 

^ number of linearly independent eigenstates $(a;) with energy E is the degeneracy of the 

energy level E. Given a Hamiltonian H, one of the fundamental questions is: What are 
^ the energy levels of the system, and what is their degeneracy? 

lO Motion of a single particle subject to an electrostatic field is usually modeled by the 

Q Hamiltonian H = — (/i^/(2m))A -|- V{x), where m is the mass of the particle, V{x) is the 

^ potential of the electrostatic field, and A is the Laplace operator. This model, however, 

^ is inconsistent with special relativity, and leads to erroneous results at high energies 

Kjjl (see [301 IM]). A better model for particles at relativistic energies was proposed by Dirac 

•'-j in |16) . On the other hand, Dirac's equations are much more complicated, and therefore 

^ their mathematical treatment is often problematic (see |lZ]). For this reason, at least in 

a some applications (see [30j and the references therein, and [T8 | I35 | ). one is often satisfied 

with an intermediate model, with the Hamiltonian H given by 

H = Ho + V{x) = {-ffc'A + m^c^y^^ + V{x). (2) 

Here c is the speed of light, and Hq = (—h?c^A + m^c^)^/^ is often called the Klein- 
Gordon square root operator, or the quasi-relativistic Hamiltonian. We often say that H 
describes the motion of a quasi-relativistic particle in the potential V{x). The purpose 
of this article is to study the stationary states and their energies for the one- dimensional 



The authors were supported by the Polish Ministry of Science and Higher Education grant no. N N201 
373136. 

Mateusz Kwasnicki received financial support of the Foundation for Polish Science. 

1 



2 



KAMIL KALETA, MATEUSZ KWASNICKI, JACEK MALECKI 



quasi-relativistic particle in the infinite square potential well: 

Jo for X G (-a, a), 
K,eii(a;) = < . (3) 

I oo otherwise. 

In the non-relativistic setting, this model is referred to as the (one-dimensional) particle in 
the box, since the potential Kveii(a^) constrains the particle not to leave the 'box' {—a, a). 
Therefore, the problem studied in the present article can be described as the quasi- 
relativistic particle in the box. It is well-known that in this case the spectrum of H is 
purely discrete, and the energy levels form an infinite sequence Ei < E2 < E3 < ... ^ oo 
(see Preliminaries for more details). Our main result can be stated as follows. 

Theorem 1.1. All energy levels are non- degenerate, and 

„ /nn -n\ he ^ f \ \ 

E^= [ — - - ] ^ O - as 00. (4) 



V 2 8/ a V^. 

It should be emphasized that the operators Hq and H are non-local. In the non- 
relativistic case, one typically solves the equation if$ = E(^ separately in the intervals 
(—00, —a), (—a, a) and (a, 00), and then matches the boundary conditions at —a and a. 
For the non-local operator H , the equation in an interval has no sense unless the values 
of $ in the complement of the interval are known. This means that one cannot apply the 
non-relativistic method; and indeed, this interval-wise method would yield the erroneous 
asymptotic formula En = [mr /2){hc/a) -\- 0{l/n), with the phase shift — 7r/8 missing. 
Non-locality of this and similar Hamiltonians has caused considerable confusion; see [20] 
for further discussion. 

The proof of Theorem |1.1| follows closely the approach of |25t [28] . The prob- 
lem studied in [25] corresponds to a massless particle, described by the Hamiltonian 
H = hc{—d'^/dx'^y^'^ -\- Vweii{x), while |28] deals with Hamiltonians based on the frac- 
tional Laplace operator {—d'^/dx'^Y for s G (0,1). The argument relies heavily on the 
results of |27| [29] , where the spectral theory for similar Hamiltonians with a single infinite 
potential wall (that is, the potential V^sx\{x) = for x > 0, V^^wix) = 00 otherwise) was 
studied. The motivation of the articles cited in this paragraph is purely mathematical. 
For the discussion of some of their physical implications, see Preliminaries. 

We remark that the Weyl-type asymptotic formula En ~ {n7[/2)hc/a as n —t- 00 is 
known since 1959, due to a general result by R. M. Blumenthal and R. K. Getoor [3j. 
Hence, formula ^ essentially provides the second term (a phase shift of — vr/8) in the 
asymptotic expansion. In fact, a stronger result is proved, with detailed estimates of the 
error term 0{l/n). The statement of this theorem requires some additional notation. 

For n > 0, we define the phase shift function (see Figure [T]) 



IT Jo s'^-fP \ 2 2\/ l + /i2 
Then is an increasing function, with limits 7^0+ = and doo = vr/8. More precisely, 

= ^ _ M + 0(/i^) as /i ^ 0+; (6) 



TT 187r 

TT l + 21og(2/i) ^ Q /log/i 



= 8 j:^;;72 + — ] as /i ^ oo; (7) 
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Figure 1. Plot of i?^. Figure 2. Plot of (solid 

line) and n7r/2 — (mca/h)^ for 
n — 1,2,3 (dashed line) when 
mca/h — 1. For the electron 
(m = 0.511 MeV), this corre- 
sponds to the potential well of 
width 2a — 0.772 pm. In this 
case, fii « 1.295, fi2 « 2.792, 
fl3 w 4.342. 



these properties are established in Section |4j We define /i„ to be the (unique) positive 
solution of the equation (see Figure |2] for a geometric interpretation) 

"^Mn 2 

Theorem 1.2. For all n = 1,2, 

he" 
a 



mca _ 

-J- fin. 



En — mc^ J of + 1 < — ( 3mc^ H ) exp (- 

n \ a J V 



mca\ 



3h 



(9) 



Remark 1.3. Formula ^ provides much more detailed information than Thanks 
to the explicit estimate of the error term, it enables passing to a limit in various regimes. 
We discuss three examples of this kind. A formal verification is given in Section |4j 

(1) By taking m — )■ 0"^ in ([o]), one recovers the results for {—(P/dx'^Y^'^, the quasi- 
relativistic Hamiltonian of a massless particle in the infinite potential well, ob- 
tained earlier in [25J. In this case we have, as in Theorem = (mi/ 2 — 
7i/8)hc/a + 0{l/n) as n — )• oo. 

(2) In the non-relativistic limit c — t- oo, we recover the classical expression 
{mr/{2a))'^H/{2m) for the kinetic energy in the (non-relativistic) particle in the 
box model. Indeed, for a fixed n > 1, we have 

E„, - mc^ = j 7^ + ( - ] as c ^ oo. (10) 



2m 



Note that the phase shift — vr/S is absent in this case, and that in the limit, the 
energies grow as n^, as opposed to linear growth in (|4]). 
(3) In (g, the error term decays exponentially fast with c and a. On the other hand, 
fin is defined in a rather implicit way. One can try to replace it with a more 
explicit expression, at the price of a slower decay rate of the error term. For 
example, we have 



E„ — mc 




nn n f h 

V 

2 V 2 mca 



h 



mca 



1 



< 12 



he 



'ID 
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(for clarity, we write for -d^ here). This estimate can be used for large energy 
levels of particles in potential wells of large width. If, for example, /i > 0, ri — )■ oo 



and a = nT[h/{2mcjj) + o(n), then (11) and series expansion immediately give 



mc 



2^, 



+ 1 nil 



+ 



n 



See (72) in Section |4] for a slightly stronger version of (11), which is suitable also 
for non-relativistic limits c — )■ oo. 

Remark 1.4. By choosing natural units, we may assume, without losing generality, that 
h = c = m = 1. This corresponds to scaling properties of the Hamiltonian if, discussed 
in detail in Preliminaries. 



Theorem 1.2 can be used to prove the following useful bounds for the energy levels. 
Proposition 1.5. For all n = 1, 2, we have 



mc 




+ 1 < En < mc 




2 mca I V \ 2 mca 

We remark that the best bounds for Er, known before were 



+ 1. 



'121 




UTX h 



2 mca 



+ 1 + mc^ \ < En< mc 




WK h 



2 mca 



+ 1, 



(13) 



proved in [H]. Note that (12) significantly improves the lower bound. For a further 
discussion, see Preliminaries. 

By a general theory, the eigenstates $„ are real- valued and vanish outside (—a, a). The 
methods of the proofs of main theorems can be used to extract some information about 
the shape of $„. In principle, the eigenstate $„ (normalized in L^(R)) corresponding to 
the energy level En is similar to either sine or cosine function (depending on the parity 
of n). Our first result is an ^^^(R) estimate. 



Proposition 1.6. Let 

fn{x) 
fn{x) 



COs{HnX)l(-a,a)ix 



(14) 



with jln defined as in Theorem L2. Then for all n 
so that 



/ri|lL2(R) < 



mc 



h 



mm 



when n is odd, 

when n is even, 
1,2, the sign of can be chosen 

10 55h\/n^, 

(15) 



n' mca 



Furthermore, is an even function when n is odd, and is an odd function when n 
is even. 



We remark that the right-hand side of (15) is 0{l/^/n) as n — )■ oo, and 0{l/^/c) as 
c — )■ oo. In particular, we recover the classical formula = fn in the non-relativistic 



limit c oo. The error in (15) comes mostly from the wrong boundary behavior of the 



approximation fn- By a suitable modification of fn near ±a, the decay rate of the error 



term in (15) can be significantly improved; see Proposition 4.15 in Section 4 for details. 



The next result gives a uniform pointwise bound for the eigenfunctions. Noteworthy, 
the results for a single potential wall in [27t .29] suggest that the best admissible constant 
in (16) is likely to be much closer to 1. 
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Proposition 1.7. For all n = 1, 2, 

ll*n|Loc(R) < (16) 



a 

The remaining part of tlie article is divided into three sections. The article is meant 
to be accessible both for readers interested in mathematical physics and those working 
in theoretical mathematics. For this reason, in Preliminaries, we discuss the formal 
statement of the problem, its various reformulations, and relation to the theory of Markov 
semigroups. Furthermore, we recall known results used later in the proofs. In Section |3] 
we present the results of f27\ 122] and discuss their physical interpretation. Finally, all 
results are proved in Section |4j 

The proofs of main theorems are rather technical. For this reason, we conclude the 
introduction with an intuitive description of the proof of Theorem In [27J, the 

(generalized, or continuum) eigenstates F^lx) were found for the quasi-relativistic particle 
with an infinite potential wall, V„i^ii{x) = for a; > and V^g^ilx) = oo otherwise; see 
Section [3] for more details. The case of the potential well, studied in this article, can be 
thought as a combination of two potential walls, one at —a, and the other one at a. We 
use the eigenfunctions F^{a + x) and -F^(a — x) corresponding to these single potential 
walls, to define an approximate eigenstate for the case of the potential well. (This 



approximation is different from the function /„ in Proposition [1.6 ) 

More precisely, we define a function in such a way that $n(a;) = a~^^'^F^(a + x) for 
X close to —a, and $„(x) = (— l)"~^a~^/^F^(a — x) for x close to a. Inside (—a, a), and 
away from the boundary, $„(x) behaves as a sine or a cosine function, in a similar way as 



in Propositio n [1.6 We prove that for a certain energy En, we have H^n{x) ~ i?ji$„(x). 
Expansion of ^n{x) in the complete orthogonal set of (true, not approximate) eigenstates 
$j(x) shows that En must be close to some (true) energy level Ek(n)- 

Careful estimates show that for each n > 3, the corresponding -E'fc(n) are distinct. 
Furthermore, by an appropriate trace estimate, it follows that there can be no other 
energy levels than Ei, E2 and -Efc(n) for n > 3, and therefore k{n) = n. Theorems 



1.1 



and 1.2 follow. Other results follow either directly from Theorem 1.2 or from intermediate 



steps of its proof. 

2. Preliminaries 

We begin with some notational remarks. Quantum states are in general complex- 
valued functions. In our case, however, all states will be real-valued. For this reason, 
unless explicitly stated otherwise, from now on all functions are assumed to take values 
in R. By a function we always mean a Borel measurable function, and we extend the 
definition of any function / defined on a subset D of R to entire R by letting f{x) =0 
for X ^ D. 

As usual, L'P{D) {p G [l,oo), D C R) denotes the space of functions f on D such 
that ||/||lp(_d) = (/f) \ f{x)\Pdxy^^ is finite, and / G L°°{D) if and only if the essential 
supremum of |/(a;)| over x G -D, denoted by ||/||loo(d), is finite. By C^{D) we denote 
the space of smooth, compactly supported functions vanishing outside D, and Co(-D) is 
the space of continuous functions in D which vanish continuously at the boundary of D 
and at ±00 (if D is unbounded). 

A function f{x) on (0,oo) is said to be completely monotone (on (0, 00)) if 
(— l)"/*-"''(x) > for all X > and n = 0,1,2,... By Bernstein's theorem (see |33| 
Theorem 1.4]), / is completely monotone if and only if it is the Laplace transform of a 
(possibly infinite) Radon measure on [0, 00). 
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We denote the Fourier transform of / G L'^(R.) by J-"/; when / G L^(R.) fl L^(R), we 
have = f(x)e~^^^dx. The inverse Fourier transform is denoted by J-"^^. Occa- 

sionally, we denote the Laplace transform of a function / by Cf, Cf{^) = f{x)e~^^dx. 
We use mostly x, y, z to denote spatial variables, and ^ for a 'Fourier space' variable. 

As usual, we write f{n) = 0{g{n)) if lim sup„_j.oo \f{n)/g{n)\ < oo, and f{n) = o{g{n)) 
when lim„^oo \f{n)/g{n)\ = Q. 

2.1. Formal definition of the Hamiltonian. In this subsection we give a formal mean- 
ing to the definition ^ of the Hamiltonian H. Further explanation of the notation used 
in formula ^ is given in Subsection 2.4 using a Feynman-Kac path integral approach. 



The free one-dimensional quasi-relativistic Hamiltonian, or the Klein-Gordon square- 
root operator, is denoted by Hq. 

2 \ 1/2 



^0= (-^V-^ + mV) . (17) 



dx 

The operator Hq is an unbounded, non-local, self-adjoint operator, acting on its domain 
T^{Hq) C L^(R), and it is defined via the Fourier transform: / G V{Hq) if and only if 
(1 -|- i'^Y^'^J^f{C) is square integrable, and 



^ifo/(0 = ynW+^^/(0- (18) 

In other words, Hq is a Fourier multiplier with symbol [h^c^^'^ + m^c^)-^/^. Note that 
V{Hq) coincides with the Sobolev space H^(R). By Plancherel's theorem, Hq is positive- 
definite, and the quadratic form associated to Hq satisfies 

{f,Hof)L^n) = ^ {J'fiO, v/^^W+^-^/(e))L^(R) > mc^ UWlHk) (19) 

for all / G T>{Hq). Formula ( [l9| ) states that the total energy, measured by Hq^ is greater 
than the rest energy mc^. 

Let V^e\\{x) denote the infinite square potential well of width 2a, as in ([S]). For brevity, 
we denote D = {—a, a); hence V^eii{x) = for x G -D and V^eiiix) = oo otherwise. 
Since V^eii{x) is highly singular, the definition of H = Hq + V^eii{x) is problematic: 
multiplication by V„c\i{x) is not well-defined on V{Hq). Clearly, Hf can only be defined 
when / vanishes almost everywhere outside D. Furthermore, when / G T>{Hq) and 
f{x) = for almost all x ^ D, then we expect that Hf{x) = H^flx) for x E D. This 
motivates the following definition. 

Definition 2.1. Let H be the unbounded, positive-definite symmetric operator with 
domain T>{H) = C^{D), D = {—a, a), given by the formula Hf{x) = HQf{x)lD{x). 
We define H to be the Friedrichs extension (or the minimal self-adjoint extension) of H. 
This means that the quadratic form of H is the closure of the quadratic form of H. More 
precisely, a function / G L?'{D) belongs to the domain T){H^^'^) of the quadratic form of 
H if and only if: 

(a) there is a sequence fk G C^{D) convergent to / in L^{D), and such that 
(/fc - //, H{fk - fi))L2(D) as /c, / oo. 

The operator H is associated to the form defined above; that is, / G L'^{D) is in the 
domain 'D{H) of H if and only if it satisfies (a), and furthermore: 

(b) there exists Hf G L'^{D) such that {Hg,f)L^D) = {g,Hf)L2(D) for g G C^{D). 
Condition (b) automatically defines the action of H on V{H). Note that in condition (a), 
we have {fk - fi,H{fk - fi))LHD) = \\Ho^^ifk - fi)\\L'HR) = \\H^^'^{fk - fdWi^D)- 
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More generally, suppose that D is an open subset of R, V{x) = for x E D and 
V{x) = oo for X ^ D. Then by Hq + V{x) we mean the Friedrichs extension in L'^{D) of 
the operator with domain C^{D), which maps a function / G C^{D) to HQf{x)lr){x). 

The above definition mimics the well-recognized non-relativistic case. There the Hamil- 
tonian corresponding to the particle in the box model is based on the Dirichlet Laplace 
operator in D, which is the Friedrichs extension (in L'^{D)) of the Laplace operator 
restricted to C^{D). Noteworthy, the self-adjoint extension of the Laplace operator 
restricted to C^{D) is not unique (an example of a different extension is the Neumann 
Laplace operator in D), hence t he ch oice of the Friedrichs extension might seem arbitrary. 

is 



On the other hand. Proposition 2.3 below states that the operator H in Definition 



2.1 



essentially self-adjoint, and hence H is the unique self-adjoint extension of H. That is 



the phrase 'Friedrichs extension of if' in Definition 2A_ can be replaced with 'the closure 
of H\ 

The following characterization of the domain of the quadratic form of H seems to be 
well-known to specialists. 

Proposition 2.2. Suppose that f G L'^{D), D = {—a, a). Then the condition (a) of 



Definition 2.1\ is equivalent to each of the following: 

(al) / G V{H^I^); 
(a2) / G V{Hl''); 

(a3) (n^c^^^ ^m^c^) 1/2 1 jr/(^) 1 2 mtegrahle on R. 

In other words, V{H^^'^) is the set of functions in V{Hq^'^) vanishing outside D. 

Note that an analogous statement is true for D = (0, oo) (with a very similar proof), 
but not for arbitrary open sets D C R (a typical example of a domain D for which the 
equivalence of (al) and (a2) fails is the complement of the Cantor set). 

Proof. Equivalence of (a2) and (a3) follows right from the fact that Hq takes a diagonal 



form in the Fourier space. Also, (a) is equivalent to (al) by Definition 2.1 and the 
definition of the square root of an operator. Finally, (a) clearly implies (a2). Hence, it 
suffices to prove that (a3) implies (a). 

For an arbitrary function / G L'^{D), denote Mirf{x) = /((I + e)^^x). Suppose that 
(a3) is satisfied. Fix h G C^(R) satisfying h{x) = for x ^ {—a, a), h{x) > for 
X G (—a, a), and J h{x)dx = 1. Let /i^ be the approximation to the identity, given by 
hs{x) = e~^h{e~^x). Note that is supported in {—ea,ea). We define = * f and 
/e = Mi^g^. Since ge G C^(R) is supported in a compact subset of (—(1 + £:)a, (1 + e)a), 
we have G C^{D). By Plancherel's theorem, 

2 1 

Wfe - /lli.(R) + - /) = TT- / - ^fm'd^, (20) 

L^(R) zvr 

where v{^) = 1 + {h^c^^'^ + m^c'*)^/^. We claim that the right-hand side converges to 
as £ —7- 0^. Once this is proved, the sequence fijk converges to / in L'^{D) and 

(/i/fc - fi/h H{fi/k - /i/i))L2(D) = \\Hl''^{fiik - /i/i)ili2(R,) as A;, / oo, that is, /i/^. 
satisfies condition (a). Hence, it remains to prove the claim. 
Observe that 



FfM) = (1 + e)FgM + 6)0 = (1 + e)Fh{e{l + e)0^f{{l + e)0. 
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and that is a continuous function. Hence J-'fe{0 converges pointwise to 
Furthermore, we have < J^h{^) < 1 and < v{^) < v{{l + s)^), so that 

< 2v{0 1 (1 + e)Th{e{l + £)0^/((l + e)0\' + 2v{0\J'f{0\' 

< 2(1 + e)Mil + ^)OI-^/((l + + 2viO\J'm)\'- 

By the assumption (a3), v{^)\J-'f{C,)\'^ is integrable. Therefore, the family of functions 
v{C,)\J^feiO ~ -^/(OP ^ (0)1)) is tight and uniformly integrable. This proves our 
claim. □ 

The properties of H for multi-dimensional domains D have been studied in [T2| [32] . 
which adopt a probabilistic definition of H; see Subsection 2A for further discussion. 

2.2. Scaling. A significant simplification of the spectral problem for H is possible thanks 
to the well-known scaling properties of the operators Hq and H. Below we show that 
the number of parameters can be reduced to one, and we may assume that, for example, 
c = h = m = 1, with the only free parameter a, half of the width of the potential well ([s]). 
This procedure is standard, see, for example, |12n32j. and in physical terms, this is simply 
the choice of natural units. 

Recall that D = {—a, a). In this paragraph only, we denote by H, Hq, $„, En and 
D the corresponding objects with h, c, m, a replaced by 1, 1, 1, a = h^^mca, respectively. 
By ^ (or ([T8|)), if /(x) = f{h-^mcx), then 

Hof{x) = mc? HqJ {h~^mcx) . 



In particular, this means that / G V{Ho) if and only if / G V{Ho). By Definition 
It follows that 



2.1 



^ , , mc - i"mcx\ , ^ n- , ^ 

^n{x) = W — $„ y-^j and En = mc^En- (21) 

We have thus proved that it suffices to consider only the one-parameter family of Hamil- 
tonians, corresponding toc = h = m = l and a > 0. 

It is easy to check that the above scaling is consistent with the results given in The- 
orems 1.1 and 1.2, and Propositions 1.5 -1.7 Therefore, it is sufficient to prove them 



only when c = h = m = 1. For easier reference, we state the simplified versions of main 
theorems. 



Theorem 1.1 '. Let En be the sequence of energy levels of the Hamiltonian 



i/=(^-^ + lJ +K.cii(x), (22) 

where KveiK^;) = for {—a, a) and V^e\\{x) = oo otherwise. Then En are non- degenerate, 
and 

En = — - — + 0{-) asn^oo. 23 

2a 8a \n I 



Theorem |1.2f . For n > 1, let fin be the unique solution of 

nvr 
~2 



^f^n = — - «/^n, (24) 
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where 




Then 

^.-V«TT|<^(3+l)exp(-^). (25) 

To simplify the notation, whenever possible, we assume that c = h = m = 1; for 
example, we prove Theorems |1. If and L2 rather than their original versions. The general. 



non-scaled case only needs to be considered in the verification of Remark 1.3 



2.3. Free Hamiltonian and its Markov semigroup. The evolution of a quantum 
state is described by the group of unitary operators exp{itH) generated by the Hamil- 
tonian H. However, when studying the spectral theory of H, it is often more convenient 
to work with the semigroup of self-adjoint contractions exp{—tH), t > (note that all 
of the operators H, exp{itH) and exp{—tH) have the same eigenfunctions) . It turns out 
that, as in the non-relativistic case, this semigroup is a sub-Markov semigroup, related 
to a certain stochastic process. Using this semigroup approach, one can relatively easily 
obtain some simple properties of H and its eigenvalues. Below we only state the results, 
and we refer the reader to [71 Chapter 5] and to [lOj for the full exposition. 

In this subsection, we focus on the free Hamiltonian (in natural units c = h = m = 1) 



1/2 



^^oH-^ + ll . (26) 



and the closely related operator 



where for brevity we let 



w{0 = Vi + ~^-^- 



Note that Aq measures the total energy (given by the Hamiltonian Hq) minus the rest 
energy mc^ = 1. For this reason, it is often called the kinetic energy operator. Both 



Hq and Aq have the same domain, and are nonnegative-definite on L'^(R,) (cf. (19)). 
Their spectra are a{Ho) = [l,oo) and (t{Ao) = [0,oo), and both operators are Fourier 
multipliers with symbols, respectively, (^^ + 1)^/^ and w{^^) = (^^ + 1)^/^ — 1. 

The operator —Aq generates a strongly continuous semigroup of self-adjoint contrac- 
tions 

To{t) = exp(-tAo), 

with t > 0; T'o(O) is the identity operator, To(t) is the Fourier multiplier with symbol 

exp(— tw(^^)), and 



V{Ao) = |/ e ^'(R) : Hm '^^^ ^ exists in L'(R) 

-Aof = hm . 

t^o+ t 



(27) 
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For t > 0, the operator To(t) is a convolution operator, and it has an everywhere positive, 
continuous, symmetric kernel function To(t; x — y): for all / G L^(R), 

/oo 
Toit;x-y)fiy)dy. (28) 
-oo 

The Fourier transform (in the x variable) of To(t;x) is exp(— tw(^^)), and To(t;x) is 
integrable function, with total mass 1. Hence, To(t) is a Markov operator, and formula ( pSj ) 
defines a contraction on every L^(R) {p G [l,oo]), and also on Co(R), the space of 
continuous functions on R convergent to at ±oo. In each of these Banach spaces, we 



can define the generator of the semigroup To(t) in a similar way as in (27); for example, 

(29) 



V{Ao;Co{R)) = if e Co(R) : lim '^^^^ ^ exists in Co{R) 



-Aof = \im . 

t-s>0+ t 

We abuse the notation and use the same symbol —Aq for the generator of the semigroup 
To{t) in any of the spaces L^(R) or Co(R). This causes no confusion: if / G 'DIAq; L^(R)) 
for different values of p, then AqJ does not depend on the choice of p. Also, if / G 
V{Ao;LP(R)) nV{Ao;Co(R)), then both definitions of Aof coincide. Finally, C^°°(R) 
is contained in V{Aq, L^CR)) {p G [l,oo]) and in V{Aq;C()(R,)), and it is the core of 
Aq in each of these Banach spaces except L°°(R). Whenever we write V{Aq), we mean 
V{Ao;L'{R)). 

For / G C^(R) we have, by Fourier inversion formula. 



Aofix) = pv / if{x) -f{x + z)Mz)dz, (30) 

J — oo 

where 

.W^i^. (31) 

TT \Z\ 

Ki is the modified Bessel function of the second kind, and 'pvj' denotes the Cauchy 
principal value integral: 



pv / {f{x) — f{x + z))i'{z)dz = lim / {f{x) — f{x + z))i'{z)dz. 

) iR\{-e,e) 



In particular, Aq is a non-local operator; see |32| Lemma 2] for a proof of (30) and (31). 

We choose, following |10) and others, to work with the semigroup generated by —Aq 
rather than —Hq, because the constant 1 is an invariant function for To(t) = exp(— tAg). 
This choice is not essential, since the relation exp(— tifg) = e~*To(t) enables one to switch 
freely between the two semigroups. 

2.4. Relativistic Levy process. Recall that the Hamiltonian for the quasi-relativistic 



particle in the box model is given by if = Hq + V^q\\{x), formally defined in Definition 2.1 
We define in a similar manner A = Aq + KveiK^^), and we let T{t) = exp{—tA). It is easy 
to see that A = H — 1 and T{t) = e* exp{—tH). We denote, as usual, D = (—a, a). 

In this subsection, we give an alternative construction of H, A and T{t), involving 
some theory of Markov processes. Probabilistic concepts described here are never used 
outside this subsection. Hence, the reader unfamiliar with Markov processes may freely 
skip this part, noting only that the operators T(t) admit a kernel function T{t; x, y) such 
that < T(t; x, y) < TQ{t; x — y) and T{t; x,y) = whenever x or ?/ is outside D. Note 
that here D = {—a, a), but the extension to more general open sets (in particular, to 
D = (0, oo), corresponding to a single potential wall Kvaii(a;)) is straightforward. 
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The operator —Aq (acting on the Co(R) space) is the Feller generator of a stochastic 
process X{t), called the relativistic Levy process, or simply relativistic process. It is a Levy 
process: it has independent and stationary increments, and cadlag (right-continuous, with 
left-limits) paths. In fact, X{t) is a subordinate Brownian motion; see [H Chapter 5]. 
We write P^^ for the probability law for the process starting at x G R, and E^. for the 
corresponding expectation. Hence, P^(X(t) G dy) = TQ{t; x — y)dy] that is, To{t]x — y) is 
the transition density of X{t). The paths oi X{t) are discontinuous: X{t) is a jump-type 



process (the intensity of jumps, or the Levy measure, is given by i'{y)dy, see (31)). Let 
D = (—a, a), and define the time of first exit from D 

td = inf {t > : X{t) ^ D} . 

Let X^{t) be the part of X{t) m D,ot X{t) killed upon leaving D: the (sub-Mar kovian) 
stochastic process with lifetime To, equal to X{t) for t G [0,Td)- Then for t > and 
X,?/ G R, the transition density of Xf^ is given by Hunt's formula, 

T{t;x,y) = To{t;x- y) - E^.(To(t - Td;X{td) -y);t> td). 

It is known that T(t;x,y) = T{t;y,x), T{t;x,y) = whenever x ^ D ot y ^ D, and 
< T{t; X, y) < To(t; x, y). The corresponding semigroup of sub-Markov operators is 

/oo 
Tit;x,y)fiy)dy = E,if{Xit))h^r^) 
-oo 

for / G L^(R) (p G [1, oo]) or / G Co(R). Since T{t)f = for any / supported outside 
D, T{t) is not strongly continuous in L^(R). It is, however, strongly continuous in L'^{D) 
(with the usual identification of L'^{D) as a subspace of L^(R)). On its L'^{D) domain, 
the generator —A of T{t) coincides with the operator —{H — 1), where H is defined 



analytically in Definition 2.1 Again, the definition of A can be extended to Lp{D) 
{p G [1, oo]) and Cq{D), with appropriate domains. We remark that T{t) maps Cq{D) to 
Cq{D) also for D = (0, oo), but not for all open sets D C R (D = R \ {0} being a simple 
counterexample) . 

Finally, the semigroup T{t) can be thought as a Feynman-Kac transform of TQ{t), 
corresponding to the potential Vweii(a:). Indeed, since {t : X (t) ^ D} Cl {td, td + e) almost 
surely has positive Lebesgue measure, the functional jl^V^c\\{X{s))ds is almost surely 
equal to zero when t < t^,, and infinity otherwise. (The last sentence remains true for 
D = (0,oo), which corresponds to the single potential wall Kvaii(2;); however, it is not 
true for general open sets D.) Hence, 

T{t)f{x) = (^exp ^ V^,n{X{s))ds^ f{X{t))^ , (32) 

with exp(— oo) understood to be 0. If V{x) is bounded, Aq + V{x) can be understood 
literally (because multiplication by V{x) is well-defined on V{Aq)), and a formula similar 
to (32) defines a semigroup generated by Aq + V{x) (see [TUt [2S])- In the case of the 



potential well, Aq + V^c\i{x) is formally not well-defined, but (32) still makes sense, and 
it is natural to understand Aq + Vy^c\i{x) as the generator of the semigroup (32). This 
explains the notation used in (22) and Definition |2.1 



We turn back to the analytical setting, and refer the interested reader to [T2l |26| [32] for 
recent developments on the relativistic process X{t), or to |22] and the references therein 
for research related to processes X{t) corresponding to more general Bernstein functions 
w{^). For various estimates of T{t;x,y) and the Green function for A in the quasi- 
relativistic setting (for multi-dimensional domains D), see [3^ 112] . Similar estimates for 
some other Hamiltonians Hq can be found in [51 El IHl ED 125] . 
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2.5. Eigenfunctions and eigenvalues. Some spectral properties of A and H follow 
immediately from the general theory of strongly continuous semigroups of compact self- 
adjoint Markov operators. There is a complete orthonormal set of eigenfunctions $„ G 
L^{D) of T{t), and the corresponding eigenvalues have the form e"*'^", where the sequence 
A„ is nondecreasing and unbounded. Furthermore, $„ G L^{D) (because the kernel 
T{t]x,y) is bounded for each t > 0), and even $„ G Co{D) (by an argument used for 
example in [15j for the non-relativistic case). By the symmetry of D = {—a, a), we have 
T{t;x,y) = T{t; —X, —y), and hence the spaces of odd and even L'^{D) functions are 
invariant under the action of T{t). Therefore, we may assume that every $„ is either 
an odd or even function. The ground state eigenvalue Ai is positive and simple, and 
the corresponding ground state eigenf unction has constant sign in we choose it to be 
positive in D. The functions are also the eigenfunctions of A, the generator of the 
semigroup T(t), and H = A + 1; the corresponding eigenvalues are A„ and En = \n + 1, 
respectively. Therefore, we have 

< Ai < A2 < A3 < ... — )■ 00, or, equivalently, 1 < Ei < E2 < E^, < ■■■ ^ 00. 

Relatively little is known about A^, En and $„, in particular, no closed-form formulae 
are available. The asymptotic result A^ ~ mi /{2a), already mentioned in Introduction, 
follows directly from a very general result of [3J (see Theorem 2.3 therein). Best known 
estimates of A„ are found in [T^ Theorem 4.4], where it is proved, also in a much more 
general setting, that: 




- 



2a) - \ \2aJ 



^ +1-1 < A„< J — +1-1. (33) 



In the non-natural units, the above estimate takes the form stated in Introduction in (13). 
The upper bound in ( [33| is a rather simple consequence of operator monotonicity of the 
function w(0 = + ly^^ — 1, while the lower bound requires much more effort. In the 
proof of Theorems |1.1| and 1.2 below we only use the upper bound of (33) for n < 2 to 



show simplicity of A2. We also apply (33) for small n in the proof of Proposition 1.5 It 



should be pointed out that for large n or large a, the estimates in Theorem 1.2 are much 



better than both (12) and (13). They also improve bounds for sums of eigenvalues 
obtained in 1191. 



2.6. Essential self-adjointness. The following interesting result seems to be new. It 
provides an important property of the quasi-relativistic particle in the box model, which 
contrasts sharply with the non-relativistic setting. For this reason we decided to include 
it, despite the fact that it is not required for the main results of the article (it is referred 



to only in the remark following Definition 2.1) 



Proposition 2.3. The operator H given in Definition 2J_ is essentially self-adjoint on 
L\D). 

Proof. Denote by Jq = {—d^/dx'^Y^'^ the one- dimensional fractional Laplace operator. 



and let J and J be defined as in the first part of Definition 2A_^ but with Hq replaced 
by Jq. On C^(R), the operator Jq — Hq i?, a. Fourier multiplier with bounded symbol 
1^1 — (^^ + 1)^^^, and hence it extends to a bounded operator on L^(R). Therefore, also 
J — H extends to a bounded operator on L'^{D). Hence J is essentially self-adjoint if 
and only if H is essentially self-adjoint. Below we prove that J is essentially self-adjoint. 
(With some effort, the following argument could be modified to work directly with H; 
however, the theory for the fractional Laplace operator Jq is much better developed than 
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the corresponding one for Hq, so references to hterature are much easier in the case of 

Clearly, J is a symmetric unbounded operator on L^{D), with domain T>{J) = C^{D). 
Furthermore, there is £ > such that || J/||l2(£)) > £:||/||l2(z)) for / G T>{J) (because the 
ground state eigenvalue of J is positive; see jll [21 |2ll [13]). Hence, is a bounded 
operator defined on the range of J. We claim that the range of J is dense in L'^{D). 
Once this is proved, it follows that the closure of is a bounded self-adjoint operator 
on entire L'^{D). This means that the closure of J is the inverse of a bounded self-adjoint 
operator, and the desired result follows. 

Suppose, contrary to our claim, that there is a nonzero g G L'^{D) orthogonal to the 
range of J. Then for any test function / G C^{D), 

= {g,Jf)L^D) = {g, ■Jof)L^R)- 

By Plancherel's theorem and Fourier inversion formula, the right term is equal to 
^o(/>^7)(0), where f{y) = f{-y) (note that f^gE C-(R) C V{Jo)). Replacing / 
by /, we obtain that for all / G C'^{D), 

Mf*gm = o. 

Let G C^{D) [e G (0,a/2)) be an approximation to the identity: supp/i^ C [— 

/le > and J hs{x)dx = 1. Furthermore, we assume that < 1/e and hs{—y) = hs{y). 

For X G (—a + e,a — e) the function f{y) = h^^x + y) is in C^{D). Hence, 

o = Mf*9m = Mhs*g){x). 

That is, Jo{he * g){x) = for x G (—a + e, a — e). 

In the remaining part of the proof we use the notion of harmonicity with respect to 
Jo (or with respect to the corresponding stochastic process, the symmetric 1-stable Levy 
process), explained in detail, for example, in [B] (see [121 [32] for the closely related notion 
of harmonicity with respect to Hq). Since in this article the concept of harmonicity is 
used only in this proof, we are very brief here. A function h is said to be Jo-harmonic 
in a domain f/ C R if for all intervals [xq — r,Xo + r] C U and all x G {xq — r,Xo + r) it 
satisfies a mean-value property 



h{x) = / h{y)P,^,,^r{x,y)dy, 

JlV\ [xo-r,xo+r] 

where Pxg^r{x,y) is the Poisson kernel for Jq. A closed-form formula for PxQ^j.[x,y) is 
known due to M. Riesz (|31|). we have 




-Pxo,r {x, y) 

Here |x — XqI < r < \y — XqI; if this condition is not satisfied, we take PxQ^r{x-,y) = 0. 

It is proved in [Gj Theorem 3.12], that weak harmonicity of hi;*g (that is, the condition 
<^(^£ * 9){.x) = for X G (—a + £, a — e)) implies (strong) harmonicity. Hence, for all 
X G (—a + 2e,a — 2e), 



he*gix)= Po,a~2eix,y)h,* g{y)dy 

J'R\(-a+2e,a-2e) 



Po,a-2eix, y)K{y - z)g{z)dzdy. 

R\(-a+2e,a-2e) J {-a, a) 
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By the exphcit formula for the Poisson kernel, 



, , 1 2a(a-2e-\x\) 1 V2 1 1 

Po,a-2e{x,y) < 



IT 



a{\y\ - {a - 2e) a - 2e - \x\ n ^a-2e- \x\ ^J\y\ - (a - 2e) 



The properties of and a short calculation yield that for some constant C, 
r CI 

/ Po^a-2e{x,y)he{y - Z)dy < = —= lR\(_a+3j,a_3j)(z); 

Jn\{-a+2s,a-2e) y a — 2e — \x\ 



see [U El E] for a similar regularization of the Poisson kernel. We conclude that 



\he*g{,x)\< I I Po,a-2eix,y)he{y - z)dy \ \g{z)\dz 

J{-a,a) \JR\(-a+2e,a-2e) / 

C If 

<^==^^^^/ \9{z)\lR.\{-a+3e,a-3e){z)dz. 

ya - 2e - \x\ J{-a,a) 

By Schwarz inequality, 

C If 

\he * 9{x)\ < , ^ = ^ / \9{z)\dz 

'a — 2s — \X\ y£ J {-a,a)\{-a+3e,a-3e) 



•J a - 2e - \X\ \J {~a,a)\(-a+3e,a-3e) 



Let e = 1/n and n — t- oo. The left-hand side converges to \g{x)\ (for almost all x G -D), 
and the right-hand side tends to 0. It follows that g{x) = for almost all x E D, contrary 
to our assumption that g is nonzero. The proposition is proved. □ 

Noteworthy, the above argument works also for the fractional Laplace operator 
{—d'^/dx'^y and the related operator {—d^/dx'^ + 1)*, but only when s G (0,1/2]. For 



s G (1/2,1], Proposition 2.3 fails: the corresponding operators are not essentially self- 
adjoint on L'^{D). 

3. QUASI-RELATIVISTIC PARTICLE IN THE PRESENCE OF A SINGLE INFINITE 

POTENTIAL WALL 

The main subject of this article is the quasi-relativistic particle in an infinite potential 
well. As described in Introduction, a potential well can be viewed as a superposition of 
two infinite potential walls, and the case of a single infinite potential wall was recently 
solved in [27112^ . The framework of these papers is purely mathematical. In this section 
we state some of the results of [27t [29] in the context of quasi-relativistic particles. 

3.1. Non-relativistic case. It seems instructive to begin with the well- understood clas- 
sical, non-relativistic model, where the free Hamiltonian is given by Hq = —d'^/dx'^ (for 
simplicity, we omit the multiplicative constant h/{2m)). The eigenfunctions of Hq are 
the free wave solutions e^^^ G R), corresponding to energies E = fi'^. These eigen- 
functions are not L'^(R,) functions. Nevertheless, any (complex- valued) ^^^(R) function $ 
can be written as a 'mixture' of free waves (or a 'wave packet') by means of the Fourier 
transform, at least when is integrable: 

$(x) = — / J^$(/i)e^'^^d/i. 
27r J-oo 

On its domain, Hq acts under the integral sign, 

1 

Ho<l>{x) = — / n^J^<l>{fi)e'^''dfi, 
27r 
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at least when n'^J-'^{fi) is integrable. In other words, the Fourier transform diagonalizes 
the action of Hq: for any $ G V{Hq), 

^(i/o*)(/x) = /x2^$(/i). 

For real-valued functions, one can replace complex-valued free waves e*^^ (/i G R) by 
real- valued ones, sin(/ix) and cos(/ix) (/^ > 0). 

Consider now a non-relativistic model with an infinite potential wall Kvaii(a;) = for 
X > 0, Vwaii(2^) = oo otherwise. The corresponding Hamiltonian -ff„aii = Hq + V„aii(a;) is 



defined as in the second part of Definition 2.1 (with non-relativistic Hq and D = (0, oo)) 



and it is the Dirichlet Laplace operator in (0, oo). The eigenfunctions of -ff„aii must 
vanish in (— oo, 0] and they are eigenfunctions of Hq on (0, oo) (here we use the locality 
of the operator Hq). Therefore, F^{x) = sin(/ix)l(o,c)o)(a^) (a* > 0) is the eigenfunction of 
-f^waii corresponding to the energy E = [j?. Again, ^ -^^((0, oo)), but the Fourier sine 
transform 

/■oo roo 

^sin/(yu) = / f{x)sin{^x)dx = / f{x)F^{x)dx 
Jo Jo 

diagonalizes the action of i^waii- More precisely, (2/7r)^/^J-'sin extends to a unitary operator 
on L2((0, oo)), and for all $ G V{H), 

J^sm(i^wall$)(/U) = /i^J'sin$(yU)- 

The free wave solutions both for Hq and for i^waii sire not square-integrable, so that, 
strictly speaking, they are not eigenstates of Hq and -ff„aii- Instead, they are generalized 
eigenfunctions, or continuum eigenfunctions, of Hq and -/Ywaii, and the integral decompo- 
sition of any state into the free waves is called generalized eigenfunction expansion. The 
primary goal of [271 129] was to provide similar generalized eigenfunction expansion for a 
class of Hamiltonians -ffwaii = Hq + V„aii(a^) with a single infinite potential wall, based on 
non-local free Hamiltonians Hq. 

3.2. Quasi-relativistic case. We return to the quasi-relativistic setting. We assume 
that c = h = m = 1 (scaling), so that Hq = (—d'^/dx'^ + 1)^/^ is a Fourier multiplier 
with symbol (^^ + 1)^/^. As in the non-relativistic case, the free wave solutions e*^^ for 
/i G R (or cos(/ix) and sin(yux) for /i > 0, when one prefers real- valued functions), are 
the generalized eigenfunctions of Hq. The only difference here is in the corresponding 
energies, which are now given hy E = (/i^ + 1)^/^ instead of the classical expression 
E = 

With the presence of the infinite potential wall Kvaii(a;) = for a; > 0, V„s.ii{x) = oo 
otherwise, the situation is more complicated due to non- locality of Hq. Nevertheless, it 
can be expected that the effect of non-locality decays as the distance to the potential 
wall increases, that is, the generalized eigenfunction corresponding to the energy E = 
[fi^ + iy^'^ should be asymptotically equal to sin(/ix+'(9^) as x — )■ oo, for some -(9^ G [0, 2tv). 
This is indeed the case, and the phase shift is given by ([S]). 



The operator ifwaii = Hq + V„i^ii{x) is defined in the second part of Definition 2.1 , with 



D = (0, oo). The results of Subsection 2.4 extend to the case of a single potential wall; in 
particular, i?„aii acts (as an unbounded operator) on any of the spaces L^((0,oo)) {p G 
[l,oo]) or Co((0,oo)), with the corresponding domain denoted by ©(ifwaii; -^^((0, oo))) 



etc. We remark that also Proposition 2^ holds true for the single potential wall, but the 
proof requires some modifications (roughly, one needs to consider Jq + 1 instead of Jq). 

As in Preliminaries, we define the kinetic energy operators Aq = Hq — 1 and A^aii = 
Aq + Vwaii(^) = -f^waii ~ 1- The following definition and theorem is a compilation of [27\ 
Theorems 1.1, 1.3 and Example 6.2] and [2^ Theorem 1.10], applied to the operator 
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Figure 3. Plot of Ff,{x/n) 
for fj. = 1/20,1/2,1,10 (blue, 
magenta, yellow, green). 



Figure 4. Plot of G^{x/fi) 
for n = 1/20,1/2,1,10 (blue, 
magenta, yellow, green). 



monotone function w{^) = (^+ 1)^/^ - 1. Note that in [23 EH], the funct ion w is denoted 
by the symbol tp. 

Definition 3.1. For fx > 0, let 

F^{x) = sin(/ix + t9^)l(o,oo)(a;) - Gf,{x), (34) 

where is given by ([s]) and is a completely monotone function on (0, oo), given by 
the formula 



Gf,{x) = [ I e '''^-ff,{r)dr ) l(o,oo)(a;), 

with 



^f^i^) = ^ /— ^ „2 I ^2 / 7^rrT2 (35) 



Theorem 3.2. The functions are generalized eigenf unctions of the Hamiltonian H^^ii 
and the kinetic energy operator A„aii for the quasi- relativistic particle in the presence 
of a single infinite potential wall; the corresponding energies are E = {(j? + 1)-^/^ (total 
energy) and A = {fj? + 1)^/^ — 1 (kinetic energy). The operators i^waii md A^aii admit the 
generalized eigenfunction expansion. More precisely, let 



POO 

n/(/i) = / fix)F^ix)dx 
Jo 



for f E L^i^) n L^(R). Then (2/7r)^/^n extends to a unitary mapping on L^((0, oo)). 
Furthermore, for $ G -^^^((0, oo)), 

<^eV{H^,u) ^ V/i2 + in$(/i) eL2((o,cx))), 

and if ^ E V{H„aii), then 

n(i7waii*)(/x) = v//i2 + in$(/i), n(A„aii$) = (V/i^ + 1-1) n$(/i). 

Furthermore, F^ E I?(if„aii, ^°°((0, oo))) /or a// /i > 0, and i/waiii^^ = (/U^ + 1)^^^^^/., 

v4walli^^ = ((/i2 + 1)1/2 



The derivation of the formula for F^ relies on a Wiener-Hopf method, and the proof of 
Theorem 3.2 involves complex variable methods; see ^27\ for details. 
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As a part of |27t Theorem 1.1] we have G [0,7r/2). By pTf Proposition 4.17], is 
continuous and different iable in > 0. The formula for the Laplace transform of F^{x) 
is also given in [27, Theorem 1.1], for ^ G C with Re^ > 0, 



F„{x)e ^^dx = r.^ ^r. exp 



2 + e 



TT 



eiog (l + 

— v+y/^ ^H- (36) 



We need an estimate of CF^{^). For real ^ this is proved in [29, Corollary 5.1], which in 
our case gives 



CF.iO < C 



i+e 
i+/i2 



(37) 



for ^ > (here we can take C = V2; a reverse inequality also holds true with a dif- 
ferent constant, but we only need the upper bound.) For complex ^ with Re^ > 0, by 
combining |271 Proposition 2.21(c) and Remark 4.12], we obtain 

\f^'+e\mo\<cifi'+\^\')CF,m) 

(with constant C = a/2). Here we use the fact that (/i^ + ^^)£F^(^) is a complete 



Bernstein function, see [27, Proposition 2.19 and Lemma 3.8]. This and (37) give 



when ^ G C, Re^ > (we can take C = 2 here). 



1±M 



(38) 



4. Proofs 

This section is modeled after [251 Appendix C and Sections 8-10]. Some of the results 
of this part do not depend on detailed properties of Hq, and can be applied to a wider 
class of models based on the free Hamiltonians Hq = w{—(P/dx'^) + C, where w is an 
operator monotone function (also called complete Bernstein function) and C G R. This 
more general case will be studied in a separate article. 

4.1. Estimates for ^o- Recall that Aq = Hq-1 = {-d"^ /dx^ + iy/"^ , and for / G C^(R), 



by (30) we have 



Aofix) = pv / (/(x) - f{y))v{x - y)dy 



{2f{x)-f{x + z)-f{x~z))v{z)dz, 



(39) 



where vi^z) is a symmetric, unimodal, positive kernel function given by (31). We de- 



note the right-hand side by ^0/(2;) (with a calligraphic letter A) whenever the integral 
converges. 

Proposition 4.1. Let x G R, 6 > 0, and let g have an absolutely continuous derivative 
in {x — b,x + b) . Then 



\Aogix)\ < I sup \g"iy)\ 

ye{x-b,x+b) 



z'^v{z)dz+ I {\g{x)\ + \g{y)\)v{y - x)dy. 

J'El\{x-b,x+b) 



18 



KAMIL KALETA, MATEUSZ KWASNICKI, JACEK MALECKI 



□ 



q{x) 



(40) 



Proof. The result follows directly from the definition of Aq (see (39)) and Taylor's expan 
sion for g, 

I'^gix) - g{x + z) - g{x - z)\ < z"^ sup \g"{y) \ 

y(^{x—b,x+b) 

when z G (—6, b). 

For 6 > 0, we define, as in [251 Appendix C], an auxiliary function: 

for X G {—oo, —b], 
(l/2)(a;/6+ 1)2 for xG [-6,0], 

1 - (l/2)(x/6- 1)2 for xG [0,6], 
1 for X G [b, oo). 

Note that q is C^, q' is absolutely continuous, < q"{x) < 1/6^ (x G R \ {—6,0,6}), and 
g(x) + q{—x) = 1. 

Proposition 4.2. Let b > 0, let f ^ L'^(R), and suppose that f"{x) exists and is contin- 
uous for X G [—6, 6] . Define 

POO 

M_i= / |/(x)|dx, Mo= sup |/(x)|, Mi= sup |/'(x)|, = sup |r(x)|, 

^0 x&[-b,b] x&[-b,b] x&[-b,b] 

Let g(x) be given by (40), and define g{x) = g(x)/(x). For x G (— oo,0), we have 

\Aog{x)\ < C{b, z/)(M_i + Mo + Ml + M^). 
More precisely, for x G (— oo, —6] we have 

Mn '•^^ 



\Aogix)\ < 



262 ^0 



and for x G (—6, 0), 

/Mo 2 Ml 



\Aogix)\ < 



V 62 



+ M2 



22z/(z)rf2 + z/(26)M_i 



^2/7(^)^2 + 2Mo / ij{z)dz + ij{b)M^i. 



Proof. We have 5f(?/) = q{y)f{y) = for ?/ < —6. Furthermore, g(?/), and hence also 
g{y), have absolutely continuous derivative in [—6,6]. For y G (—00,6) \ {—6,0} we have 

\q{y)\ < 1, \q'{y)\ < 1/6, \q"iy)\ < 1/62, ^.(^) ^ /(?/)g"(?/) + 2/'(?/)g'(y) + r(y)g(i/). 

It follows that for y G (—00,6) \ {—6,0}, 

, Mo 2Mi 

\9"{y)\ <^ + -^ + M,. 

Let X G (—6,0). By Proposition 4.1, we have 
'Mo 2Mi 



\Aog{x)\ < 



b"^ ^ b 



+ M2 



z v{z)dz + 



{M^ + \g{y)\)v{y-x)dy 



i'R\{x-b,x+b) 

Since x — 6 < —6, we have (7(1/) = for y G (—00, x — 6]. Furthermore, for y G [x + 6, 00) 
we have \g{y)\ < \ f{y)\ and i/(?/ — x) < z^(6). It follows that 

'Mo 2Mi ^ 
— - + 
62 ^ 

When X G (-00,-6), then 5f(x) = and g"{x) = 0, and we simply have (again using 
monotonicity of u) 

-b 



\Aog{x)\ < 



z^u{z)dz + 2Mo / u{z)dz + u{b)M^i. 



Aog{x)\< \g{y)\u{y - x)dy < Mq I q{y)u{y - x)dy + u{2b) I \f{y)\dy. 
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Using the inequalities q{y) < {b + y)^/{2b^) and u{y — x) < v{h + y) for tlie first integral 
on the right-hand side, we obtain 

M f^^ 

\Aog{x)\<^j z^u{z)dz + u{2b)M_,. □ 

4.2. Approximate eigenfunctions. Definition and Theorem 3^ yield generalized 
eigenfunctions for the operators if^aii = Hq + V^siiix) and A^aii = + V„i,ii{x) = 

-f^waii ~ I5 corresponding to the eigenvalues (yU^ + 1)^/^ and (/i^ + 1)^/^ — 1, respectively. 
Below we construct approximations to eigenvalues and eigenfunctions of H = Hq+V„c1i{x) 
and A = Ao + V^eii{x) = H - 1. 

Recall that D = (—a, a) for some a > 0, the half of the width of the potential well ([s]). 
We fix an auxiliary number b G (0, a); later we choose 6 = a/3 to optimize constants. For 
n > 1 let /i„ be a solution of 

afin + ^ii„ = Y (41) 
(this is a scaled version of ([s])), and let 

E„ = V/i^ + l, A„ = V/i^ + 1 - 1. (42) 

Here -(9^ is defined by ([s]). Since a/i + 'd^ is a continuous function of /i, it takes all 
values between 'd^^ = and 00. In particular, a number /i„ satisfying (41) always exists. 
Furthermore, since "i?^ is nondecreasing in /i > (see the next section), the solution /i„ 
of (41 ) is unique. 

Our goal is to show that En is close to some eigenvalue of H, or, equivalently, that A„ 
is close to some eigenvalue of A. To this end, we define an approximate eigenfunction (fn 
of H and A, using the eigenfunctions Ffi^{a + x), Fjx„{a — x) for single potential walls at 
—a and a. We let (see Figure [s]) 



^n{x) = q{-x)Ff,„{a + x) - (-l)"g(x)F^„(a - x), (43) 



20 



KAMIL KALETA, MATEUSZ KWASNICKI, JACEK MALECKI 



with the auxihary function q defined by (40). Here x G R, but we have V9„(x) = for 
X ^ D, so that (fn € The notation introduced above is kept throughout this 

section. Note that A„ depends on a and n, while <^ra(x) depends also on b. Note also that 
(pn is not normed in L^(R), its norm is approximately equal to ^/a (see Lemma 4.5); this 
is why we denote it with a lowercase letter, and reserve for 'pn{x)/\\'pn\\L^{D)- 

Lemma 4.3. We have (pn € '^'(A) and A(pn{x) = Aoipn{x) for almost all x E D. 

Proof. For brevity, in this proof we write ft = ftn and if = (p^. The operator A = 
Aq + Kveii(a;) is defined in Definition 2.1 Hence, we need to check two conditions: (a) 
there are fk G C^{D) such that in L'^{D) and {fk - fi, ^o(/fc - fi))L^D) as 

/c,/ oo; and (b) {ip,Aog)L2(^D) = {Ao(f, g) l2(d) for g e C^{D). We first verify (b). 

Note that Aoip{x) is well-defined for all x G D\{—b, b}, since (f is smooth in D\{—b, b} 
and bounded on R. Furthermore, by the definition of (see (39)), for any g G C^{D), 



Ao(f{x)g{x)dx — / (f{x)Aog{x)dx 



{g{x + z)ip{x) + g{x — z)ip{x) — g{x)(p{x + z) — g{x)(p{x — z))u{z)dzdx, 

which is by Fubini, provided that the double integral converges. Denote the integrand 
by /(x, z)z/(z), and let e = dist(supp (?, R\D)/3, so that supp (7 C (— a+Se, a — 3e). When 
z > e, then \I{x, z)\ < Ci{(f, g). Suppose that z G (0, e). If x ^ {—a + 2e, a — 2e), then 
I{x,z) = 0. Otherwise, by first-order Taylor's expansion of I{x,z) around z = (note 
that I{x, 0) = {d/dz)I{x, 0) = 0) with the remainder in the integral form, we obtain that 
|/(a;, z)| < C2{'f,g)z'^. We conclude that \I{x, z)h'{z)\ < C^^ip, g)mm{l, z^)h'{z), which 
implies joint integrability of I{x, z)i'{z). Hence, condition (b) is satisfied. 

By Proposition 2.2, in order to verify (a), one only needs to show that {l+^'^Y^'^\J^(p{C)\^ 
is integrable on R. Let f{x) = q{a — x)-F^(x), so that <^(x) = f{a + x) — (— l)"'/(a — x) 
(see (43)). Hence, it suffices to prove integrability of (1 + ^'^Y^^lJ^fiO?. 

Fix e > and let q{x) = q{a — x)e'^^ . It is easy to check that the distributional 
derivatives g, q' and q" are integrable functions, and the third distributional derivative of 
q{x) is a finite signed measure on R. Hence, q{x) has the same property. Therefore, J^q{^) 
and J^q"'{i) = —ii^J^qH) are bounded functions, and so \ J^q{i)\ < C4(e,a,6)/(1 + \^\)^. 

The Fourier transform of e~^^F^(x) is CFfi{e + i$,), and the Fourier transform of f{x) = 
q{a — x)F^{x) = g(x)e~^^'F^(x) is given by the convolution 



1 

27r 



Suppose that ,^ > 0. To estimate |J-'/(^)|, we write 

1 r°° 1 

J^q{^ - s)CF^{e + is)ds + 



J'fiO 



2n 



?/2 



27r 



J^q{s)CFf,{e + i{^- s))ds. (44) 



«/2 



By (38), we have 



;i + |s|)3/2- 



Hence, 



mi 



s)CF~^{e 



< 



is)ds 
Ce{e,il)C4{e,a,b) 



< 



:i+e/2)3/2 
1 



(l + e/2)3/^ 7^/2 (l + (e-s)) 



\J^q{^-s)\ds 



C/2 

ds < 



C(i{e,^)Ci{e,a,b) 
(1 + e/2)3/2 
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For the other integral in (44), we use + is)\ < CQ{e,jl): 

J^q{s)CFf,{e + i{^- s))ds 



?/2 



Therefore, for ^ > 0, 



1-^/(61 <C7(£,a,6,/i) 



:i + ie/2|)3/2- 



Since J-'f{—^) = J^fiO^ ^^e above estimate holds for all G R. We conclude that for all 



i+eY^'\j'fm''<mie,aAf^)) 



r-.\\2 



and the right-hand side is integrable. The proof is complete. 

For brevity, in the remaining part of the article we use the following notation: 



□ 



(45) 



G^{x)dx, G^,h{.x) = G^{x — b) + G^{x + h). 



The following results are stated for A and A„; their reformulation for H and En is straight- 
forward. 

Lemma 4.4. We have 
More precisely, we have 



\\A(pn - K^n\\\2(^D) < 2(a - h) 



262 



+ 26 



(46) 



\ G-,^M - 2hG'-^^^,{a) + h^Gl^^,{a))u,{h) 

+ 2G/i„,fe(a)z/oo(6) + z^(6)//in + ^ + ~1~ 

Proof. The argument below follows exactly the method developed in |!25], Lemma 1], and 
applied also in [28^, Lemma 1]. Since in this proof we consider n fixed, we simplify the 
notation and write fx = fin, \ = \n = (/i^ + lY^'^ — I and (p = (pn- 

Recall that Fp{x) = sin(/ix + -(9^)1(0 oo)(x) — Gjx{x), and G^{x) = for x < 0. Note that 
by (HI), 

sin(/i(a + x) + -d^) = sin(^ + fix) 

= -(-1)" sin(f - fix) = -(-1)" sin(/i(a - x) + ?9^). 



Hence, by (43) and (34), for all a; G R, 

i^(x) = q{-x)F^{a + x) - (-l)"g(x)F^(a - x) 

= sm{f + flx)lD{x) - q{-x)Gfi{a + x) + (-l)"g(x)G'^(a - x) 



(47) 
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In a similar way, for all x G R, 

F^(a + x) + (-l)"F^(a - x) = sm{f + flx)lB.\D{x) sign a; 

-(G^(a + x) + (-irG^(a-x)), 



so that, again by (43), for all x G R, and 

^(x) - F^ia + x) = {qi-x) ~ l)F^(a + x) - (-l)"g(x)F^(a - x) 
= -g(x)(F^(a + x) + (-l)"F^(a - x)) 

= g(x)(G^(a + x) + (-l)''G^(a - x)) - sin(^ + /ix)l[a,oo)(a;). 

Denote 

fix) = Gi,{a + x) + (-l)"G^(a - x), 

^(x) = g(x)/(x), 

h{x) = sm{f + /ix)l[a,oo)(a:). 

It follows that (f{x) = Ffi{a + x) + g{x) — h{x). For x G (—a, 0), we have AoFfi{a + x) — 
AF^(a + x) = and h{x) = 0. Hence, for x G (—a, 0), 

\Ao^{x) - A(^(x)| < \Aog{x)\ + \Aoh{x)\ + |A(?(x)|. (48) 

We estimate each of the summands on the right-hand side separately. 

First, we apply Proposition 4.2 for the function /. Since Gfi_{x) is a completely mono- 
tone function, Gjx{x), — G^(x) and G|^(x) are nonnegative convex functions of x G (0, oo). 
Hence, 

Mo= sup \f{x)\<Gf,{a-h)+G'^{a + h) = G'^^b{a), 

xe{-b,b) 

M,= sup |/'(x)|<-G'^(a-6)-G'-(a + 6) = -G;,», 



M,= sup \f"{x)\<Gl{a-b) + Gl{a + b)=G;^,{a). 



xe(-b,i 



Furthermore, 



M_i 



\f{x)\dx < / Gp_{a + x)dx+ / Gji{a — x)dx 



G-^{y)dy = F^. 



By Proposition 4.2, for x G (—a, — 



and for x G (—6, 0), 

1^05' W I < ^— p h 2G/i,6(a)z/oo(0) -h z/(6)J^. 

Note that |5'(x)| = for x G (—a, — &), and for x G (—6,0) we have 

|A,(.)| < M^M < ^ < 



2 - 2 - 2 

Finally, z/(y) decreases for ?/ > 0. Hence, integrating by parts, we obtain that for x < 0, 

2u{a) 



\AqKx) 



sin(^ + fix)iy{y — x)dy 



< 
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By combining the above estimates and (48), we conclude that for x G {—a, —6), 



\Aq({>{x) - X(fn{x)\ < h z^(26)/^ + 
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and for x G (—6, 0), 



\Aoip{x) - V„(x)| < 



{G~,,,{a)-2hGL^,{a) + h^Gl,{a))u,{h) 



2z/(a) 



By symmetry, similar estimates hold also for x G (0,a). Integration in x gives for- 
mula (46) with Aip{x) in the left-hand side replaced by li:){x)Ao(p{x). The result follows 

□ 



by Lemma [4.3[ 
Lemma 4.5. We /iai^e 



More precisely, 



sin('(9^ 



4//i„ < ||<^n||L2(^) <a + 



(49) 



Proof. As in the previous proof, we write jl = fXn, A = A„ and (p = (pn- By (47), 

(^(x) = sin(^ + jxx)li:i{x) - q{-x)Gf,{a + x) + {-l)"-q{x)G^{a - x), 
It follows that 



I^IIl2(d) -« 



J (^(sin(f + 

+ / x)G'^(a + x) — (— l)"'g(x)G'^(a — x))^ (ix 

J —a 

-2 1 sm{^ + fix) {q{-x)Gfi{a + x) - {-l)''q{x)Gp,{a - x)) dx. 

J —a 



By direct integration and (41), 









j dx 




-f 






2 J- 



1,./ „~\ ./ sin(2'i9/i) sin('(9/i) 
- |sm(n7r — 2/ia) — sm(ri7r + 2/ia)| = — < 



4/i 



2/i 



Furthermore, 



sin(^ + jlx) {q{—x)Gfi{a + x) — (— l)"g(x)G'^(a — x)) dx 

< I {G-^{a + x) + G-^{a-x))dx <2I-^. 



24 KAMIL KALETA, MATEUSZ KWASNICKI, JACEK MALECKI 

Finally, using Gjiix) < G'/i(0+) = sin(^9^), 

0< / {q{-x)G j,{a + x) q{x)Gji{a-x)fdx 



<2( / {G/i{a + x)ydx+ / {Gj:,{a - x)ydx 

\J —a J —a / 

<2sm{}^^J Gfi{a + x)dx + J Gfi{a — x)dx^ <4/^sin'(9^, 

and the lemma is proved. □ 

Let cr{A) denote the spectrum of A. Recall that the spectra of A and H are purely 
discrete (see Subsection 2.5), and the eigenvalues of A and H (repeated according to 
multiplicity) are denoted by A„ and = A„ + 1, respectively. 

Proposition 4.6. We have 

dist(A„, aiA)) < \\^^r. -Xn^n\\LHD) ^^^^ 

Wn\\L'2{D) 

More precisely, let ^'^^^ and A"'^'^ be the restrictions of A to the (invariant) subspaces 
of L'^{D) consisting of even and odd functions, respectively. Then (50) holds with cr[A) 
replaced by cr(A'^™'^) when n is odd, and by a{A°'^'^) when n is even. 

Proof. Let E{dX) be the spectral measure of A. Since ipn G L^{D), we have 

M'^n - A„(^„||^2(£,) = / \\-\n\'^{<fn,E{d\)(^n) 

Ja{A) 

> (dist(A„,a(A)))2 / {^n,E{dX)^n) = (dist(A„, (t(A)))2||(^„||2 

Ja{A) 

as desired. Furthermore, since is either an even or odd function, we can replace A by 
^even ^o^d^ respectively, depending on the parity of n. □ 



4.3. Estimates of 'd^ and G^. Before we apply the results of the previous subsection, 
we need some detailed estimates. Recall that the (singular) kernel i/ of A is given by (31 ), 
i>{z) = i^i(|2;|)/(7r|2;|), where Ki is the modified Bessel function of the second kind. 

Proposition 4.7. For x G R \ {0}, we have 

— <H^)<- — ^ . 51 

Proof. Let a; > 0. By the definition of the Bessel function and a substitution sinht = s, 
we have 

/■OO /"OO 

e^iri(x) = / e-^(™'^*-^)coshtc/t= / e'^^^^-^^rfs. 
Jo Jo 

Since \/l + — 1 > max(0, s — 1), 

rl poo -| 

e^i^i(x) < / lds+ e-'^^'-^^ds = 1 + - . 







We conclude that e^z/(x) = e^Ki{x)/{7Tx) < (x + l)/(7ra;^), which is the desired upper 
bound. In a similar manner, \/l + s"^ — 1 < s, and therefore e^Ki(x) > 1. The lower 
bound follows. □ 
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Since (1 + x)e~^ < 1, we have, in particular, z/(x) < l/(7rx^). This simple bound is 
often used below. With the notation of (45), for x > we therefore have 

z/o(a;) = / z'^i'{z)dz < — min (x, / (1 + z)e~^dz] = — min(x,2). (52) 



Furthermore, 



for X > 0. 



r'OO 

i^oc^ix) = I iy{z)dz < I — 7^dz = — (53) 

TTZ"^ TTX 



Proposition 4.8. The function -d^ is increasing in fi, converges to as ^ ^ 0^ , and to 
n/S as oo. Furthermore, for > we have 



dd^ 2 I 

< < - — ^ . 54 

Ct/i TT 1 + /X^ 



maxlO,^- — ) <^9^<min(^,^). (55) 



In particular, 



8 nfi J \ TT 

Proof. By ^ and dominated convergence, for /i > we have (see |27l Proposition 4.17]) 



d^^ 11/"°° 1 



dr 



/ (56) 

- 1 1 1 + /i + y/l + /i^ _ 

~ vr /i(i + /i2) 1 - + yr+T? ' 

the integral above is evaluated using the first Euler substitution. In particular, in- 
creases with fi. Furthermore, ^ can be rewritten as (see [27t Proposition 4.16]) 



1 1 + l+P?/s^ 



1+^2 

Hence, by dominated convergence, — )■ as /i — )■ and 

lim ^ =- [ ^ log ^^^^^ ds = -[ ^ log -ds = ^] 

for the last equality, see \27i Proposition 4.15]. These properties of i)^ were already 
noticed in \27i Example 6.2]. By (56) and the inequality log(l + s) < s {s > 0), we also 
obtain that 

d^^ 1 1 2a 2 1 

— < -. — , < -, /i > 0, 

dfl TT /i(l + /i^) 1 - /i + ^1 + ;i2 TT 1 + Ai 

as desired. □ 

By series expansion of d-d^/dfi (using (56)) near and oo, one easily obtains ^ and ([T]). 

For any a > 0, a/i + t?^ is increasing in /i G (0,oo). Hence, as it was observed in 
the previous section (see the remark following (41)), the solution /i„ of the equation 
a/i + = n7r/2 (a > 0, n > 1 fixed) is unique. 
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Proposition 4.9. For a > and n > 1, we have 



a V 2 

In particular, 



1 /nvr \ _ 1 /rivr 

- ' i^nn/(2a) ) < /^n < " ' 



a V 2 



min(8n/ (vra), 1) 
4a2 + 7r2(n- 1/4)2 ■ 



< /in < 



2a 8a 



2a 



Proof. Since "i^u is increasing, we liave 



1 /nvr 
a \'2 



/nvr \ _ nvr 



(57) 
(58) 
(59) 



This gives (|58|) and the lower bound of (|57|). Furthermore, by the mean value theorem, 

2 '&mT/{2a) 



dsil) and (|58E 



1 /riTT 



a 



n-K l{2a) 



1 



(^9 



n-K I {2a) 



1 /nvr 

' Ain I < 



7ra 1 + /i2 V 2a 



vra^ 1 + /i; 



r,2 



This, combined with the inequality /!„ > {{n — l/4)7r/(2a)) (which follows from (59)) 
and (55), gives the upper bound of (57). 



□ 



We now turn to estimates of G^. 
Proposition 4.10. For /i, x > 0, 



G^{x) < min 



2/i 

TT 



0.383 



and 



f°° fa 
/ GJx)dx < min ( — 
Jo ' - \8 



fi 0.217 



/i 



(60) 



(61) 



Proof. By [271 Lemma 4.21], G^(x) < G^(0+) = sin(^9^). This and (55) give (60): 



2/i 



TT 



G^{x) < sin(^9^) < min I sin — , sin — < min I — , 0.383 . 



TT 



2fi 



IT 



Furthermore, again by [27, Lemma 4.21], with w{^) = (^ + 1)^^^ — 1 we have 

cos(7?^) 



G^{x)dx 



1 

1 



cos 



COS 



2v/rT7?(v/rT7?-i) 



(62) 



l+/i2 



2Vl + /i2 

We have l/-\/l + s > 1 — s/2 + s^/32 for s < 8. Hence, whenever < 8, 



'i + \/rT 



+ 



> 



1-^ + '^ = ! 

4 64 



When /i2 > 8, the right-hand side is negative, so the inequality between left and right- 
hand sides holds for all /i. It follows that G^{x)dx < /i~^(cos('(9^) — 1 + fJ^'^/S) < /i/8. 
This proves one part of (61). Furthermore, by (55) and concavity of cos on [0,7r/8], 



cos('i9^) < cos ( max 



TT 



TT/i 



< min cos(7r/8) + 



2sin(7r/8) 
7r/i 
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for all /i > 0. Hence, by (62), 



/i / G^{x)dx < min I cos(7r/8 



2sm(7r/8) 

TTyU 



1 - 



1 + 



2Jl + ^i^ 



(63) 



When /i < 7r/(2 tan(7r/8)), the minimum in (63) is equal to 1, and hence the right-hand 
side of (63) is clearly an increasing function of /i G (0, 7r/(2tan('7r/8))]. We claim that it 
is increasing also when /x > 7r/(2 tan(7r/8)). Indeed, 



d 



cos(7r/8) 



2sin(7r/8) 

TTfJ, 



2^1+^2 



2sin(7r/8) 



2x/2(l + fi-^f/^Jfi-^ + v/l + /i-2 



TT 



the right-hand side is increasing in > 7r/(2 tan(7r/8)) and positive for fi = 
7r/(2 tan(7r/8)). Hence, the derivative of the right-hand side of (63) is positive for 
fi > 7r/(2 tan(7r/8)), and our claim is proved. The right-hand side of (63) converges 
to cos(7r/8) - I/V2 as /i — )• 00. We conclude that 



G^{x)dx < — ( cos(7r/8) — 
A'' 



< 



0.217 



proving the other part of (61). 



□ 



Proposition 4.11. Let u{x) = (1 + x)e ^ . For fi,x > 0, we have 



Gf,{x) < 



0.226m(x) 

jji^x"^ 



-G'Ax) < 



0.458M(a;/2) 



jjfix'^ 



G-XXX'-^"^. (64) 



Proof. By ([35| and estimates /i/(l + /i^)"*^/^ < 1, (r^ — 1)^/^ < r, fj? + r'^ > /x^, we have 



By integration, with the notation u{x 

V2 



G^{x) 



7.(r)<— ^l(i,oo)(r). 

1 + x)e~^, 

V2il + x)e-^ V2 u{x) 



e ^^rdr < 



2tt /i^x^ 



2n jji^x"^ 



In a similar manner, using (1 + x + x^/2)e ^<2(-\/3 — l)e ^^/^■u(x/2) < 1.016'u(x/2) 
(which is easily proved by elementary calculus). 



-G'{x) 



^/2_ 
27r/i2 



" -XT 2, / 72 (1 + x + xV2)e-" 72 1.016m(x/2) 
e ""W^dr < — ^ — < „ \ 



and by {1 + x + x"^ /2 + x^ /6)e ^ < 4e v^/^m(x/2) < 1.176'u(x/2) (again easily proved by 
a short calculation). 



, 3^2 (1 + x + xV2 + xV6)e-^ 3^2 1.176m(x/2) 

G nix) < ^ — < n \ 



□ 



TT 



Notation u{x) = (1 + x)e ^ is kept throughout this section. Note that u{x) is a 
decreasing function of x G (0, 00). 
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4.4. Estimates for eigenvalues. We are now in position to apply Lemmas AA and 4^ 
Recall that A„ and E'^ = A„ + 1 are the eigenvalues of A and H, respectively, and 
< Ai < A2 < A3 < ... — )■ 00. The corresponding normalized eigenfunctions are denoted 
by $„. On the other hand, A„ = (/i^ + 1)^/^ - 1, E„ = {fxl + l)^/^ and <fn/\\^n\\mD) 
(see (41), (42) and (43)) are approximations to A„, En and $„. In the definition of <^„, 
we take b = a/3 (this choice optimizes some constants below). The following results are 
stated for A, A„ and A„, but they can be trivially reformulated for H, En and En- 



Lemma 4.12. For every a > 0, the eigenvalues A„ are simple. The eigenfunction $„ is 
an even function when n is odd, and $„ is an odd function when n is even. Furthermore, 
for all a > and n = 1,2, 

(l + a/3)e-"/3 / oQc;« inQi/i\i/2 

\\AiPn - K'fn\\L'2{D) < 



j2na 



3/2 



„ 9.356 10.314 
2.497+ ^ + — 



If a > and find > 1.251, then 
(1 + a/3)e-°/=^ 



An — A„ < 



, „ 9.356 10.314 
2.497+ ^ + — 



1/2 



1.251 

Hn^ 



(65) 



-1/2 



(66) 



Proof. The argument is divided into four steps. 

Step 1. We estimate each term on the right-hand side of (46) and (49), so that in the 
next step we can apply Lemmas 4.4 and 4.5, We use the notation (45). 

Recall that h = a/3. We introduce the function w{^) = (^+1)^/^ — 1, as in Preliminaries. 
Denote, in this step only, k = 1/3, so that b = na. Since (a — 6)/2 > b and a — b > b, 
by (64) we have 

G,,(a) = G,(a - b) + G,.(a + b) < 0.226 ' + + u{b) < "■'''^^"^ 



-G;,(«) ^ -G'^(a - ,) - G'^(a + b) < 0.458 " ' "'"^^ < 



Gf,fi{a) = G^{a - b) + G^{a + b) < 1.589 u{b) < 



Hence, 



,9^„ / ^ . 9x u(b) 2.746u(6) 

b^G'' Ja) < 0.636 + 2 ■ 1.739k + 8.548^2 < ^ ' 



By ([52]), z/o(6) < 6/7r and z/o(26) < 26/vr. Hence, 

G'^,b(a)z/o(26) ^ 0.636m(6) ^ 0.608m(&) 



262 



(6-^,5(0) - 26G;,b(a) + 6 'G;:,fe(a))^o(6) ^ 2.746tx(fe) ^ 2.623^(6) 



62 



By (53) we have Uooib) < l/(7r6), so that 



2G^,b(a)z/oo(&) < 



2-0.636u(6) 1.215u(6) 



< 



Since w{ij?) = (/i^ + 1)^/^ _ ]^ < -^g qI^q have 



ti^(yU^)G'^,6(a) 0.636m(6) 0.318m(6) 



< 



2/ia2 
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By (51), z/(x) < u{x) / [tix^). This and (61) give 



0.217u(6) 0.622u(6) 



< 



0.217m(6) 0.156u(6) 



2t/(a) ^ 2u{a) ^ 0.637M(a) 



Step 2. Recall that A„ = w(/i^). For brevity, denote = fin. By Lemma 



4.4 



we obtain 



\\ALpri - K'fnWhi^D) < 2(a " &) 



0.608m(6) 0.156u(6) 0.637u(6)\ 



+ 



+ 



'2.623m(6) 1.215m(6) 0.622m(6) 0.318m(6) 0.637m(6)Y 



jjLO? 



iuih)f ( „ 9.356 
< ^ , \ 'I I 2.497 



10.314 



This proves (65). Lemma 4.5, "i?^ < 7r/8 and (61) give 



IV'n|lL2(D) 



sin?9. _ 0.383 0.868 

> a z-^ - 4/^ > a r ^ 



1.251 



/i 



/i 



/i 



/i 



(67) 



By Proposition 4.6, there is an eigenvalue Xk{n) {O' priori k{n) may depend on a, though 
later we show that in fact k{n) = n) such that 

1/2 / 1 oc:i \ -1/2 



lA 



k(n) 



An < 



u{h) 



„ 9.356 10.314 
2.497 + + — 



1 



1.251 



(68) 



provided that the right-hand side is well-defined, that is, /i„a > 1.251. In the next step 



we show that this is the case when n > 3. Proposition 4.6 also states that we may 



assume that the eigenfunction has the same type of parity as yj^, that is, it is an 

even function for odd n, and an odd function for even n. 



For n > 3, we denote the right-hand side of (68) by e„, and define the intervals 

In (Aft Eyi', Afi -|- ^n)' 

We have thus proved that each of the intervals J„ {n > 3) contains an eigenvalue Xk{n)- 
Below we show that these intervals are mutually disjoint. 



Step 3. Clearly, decreases with n. Furthermore, by (58), 0/23 > 37r/2 — tt/8 > 4.319 



and a/i4 > 47r/2 — tt/S > 5.890. It follows that the right-hand side of (68) is well-defined 
for n > 3 (as claimed in the previous step), and 



^3 < 



0.2007rM(6) 



£4 < 



0.1287rM(6) 



a a 

Since w{fi^) = {jj? + 1)^/^ — 1 is a convex function of /x, by the mean value theorem, for 
/ii < /i2 we have 

/il(/i2 - /il) 



u;(/i2) - w(/ii) > 2/iiu;'(^i)(/i2 -/^i) n ^ ■ 

Hence, using also ( [58| and /i„ > 4.319/a {n > 3), we obtain that for n > 3, 

K{K+i-~^n) (0.200 + 0.128)7rM(6) 



(69) 



1 + At 



r/2 



> 



4.319 371 37ru{b) ^ 37r 



Va2 + 4.3192 8a 8c 



4.319 



8a VVa2 + 4.319^ 



-M(a/3) > 0; 
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2 4 6 8 10 



Figure 6. Plot of the intervals (aw((n7r/2 — 7r/8)^/a^) — a£„, a w((n7r/2)^/a^) + a£„) 
for n = 3,4,5,... (each of them contains aA„, an eigenvalue multiplied by a), and the 
function a'w{'K'^ / a?) (the upper bound for a\2, the second eigenvalue multiplied by a); 
here u;(C) = {S. + lf''^ ~ 1. 



the last inequality is proved as follows: 



- \ SJ V 36; V 3 36 108; 
and the right-hand is equal to 1 for s = 0, and is strictly decreasing in s > 0: 

— 1 + - + — + e""/^ = (s^ - 6s + 18) < 0. 

dsVv 3 36 108; ; 324 ^ ' 

It follows that the intervals /„ {n > 3) are mutually disjoint. (In fact we have proved that 
the larger intervals (w((r?,7r/2 — 7r/8)^/a^) — £:„,w((r;,7r/2)^/a^) are mutually disjoint, 
see Figure [6]) Since Xk{n) ^ ^n, all numbers k{n) [n > 3) are distinct. Furthermore, by 
an argument similar to the one used above, 

7r(/i3 - vr/a) 0.2007rM(6) 



ayl + vr^/a^ 



>-(l - 0.200n(a/3) ) > 0; 



TT 



indeed, for the last inequality we have 



u(s/3) ^ < 1 + - e-'^\ 1 + - < 1 + - e-'^' < - < 

^ ' ' TT -V 3/ V 9-V 3/ - e 0.200 

Hence, Afc(„) > w^n'^/a^) for n>3. On the other hand, by ([33]), Xj < w((j7r)^/(2a)^). It 
follows that k{n) > 3 for ri > 3. To show simplicity of eigenvalues, it remains to prove 
that there are no other eigenvalues than Ai (which is known to be simple), A2 (which lies 
in the interval (Ai, ^(vr^/a^))) and Xk{n) > 3; all distinct and greater than ^(vr^/a^)). 
This is done in the next step. 

Step 4- Recall that T(t;x,y) (x,?/ G R, t > 0) and To(t;x — y) {x,y & D, t > 0) are 
kernel functions of the (sub-)Markov operators exp{—tA) and exp(— tAg). The Fourier 
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transform of To(t; x) is exp{—tw{C,'^)), and < T{t; x, y) < To(t; x — y) for all x, y G -D = 
(—a, a), t > 0. For t > 0, we have the following trace estimate (cf. [21 [23]; |25| Section 9]; 
Section 5]) 



E 



j (^^e-^^\^j{x))'^dx = j T{t;x,x)dx 



< 



To{t;0)dx 



D 



2a 



IT 



-tw 



In the last step, Fourier inversion formula was used. Below this estimate is compared 
with the lower bound for the sum of exp(— Afc(„)t). 
For a > and n > 3, let 



r/2 



Clearly, 5n is decreasing and 5„ — )■ as n — )• oo. By (69) 
Hence, 



Furthermore, by the definition (41) of /i„ and monotonicity of ^9^, 

1 



TT 



< 



TT 



2a a ' '^"^ - 2a 

Using these two inequalities and monotonicity of 6n, we obtain that for any > 3, 



2a 



n=N 



n=N 



oo „ 

> ^{p-n+l - fj'n)exp{-tw{{fMn + SnY)) > j 

n=N 



-tw 



Ajv+5jv 



Let Jjv be the set of those j > I for which j ^ k{n) for all n> N. It follows that 



E 



E^-''*-E 

i=i 



< 



2a 



TT 







n=Af 

e-*"'^^ < — (/ijv + ^iv). 
vr 



Taking the limit as t — t- O"*", and using the definition (41) of fi^, we obtain that 

2a 



#Jn <-{f^N + SN)=N--{^ 
TT TT 



a5 



N) 



Note that as — oo, we have 6n 0, and "(9^^^ — )■ vr/S. Take sufficiently large, so 
that a6N < '^fij^ (a > is fixed). For this N, we have ^^Jn < N. On the other hand, Jjv 
must contain Ai, A2 and Afc(„) for n = 3, 4, — 1. It follows that # Jat = A^ — 1, and 
therefore there are no other eigenvalues than Ai, A2 and Xk{n) > 3), as claimed in the 
previous step. In particular, k{n) = n for all n > 3, and all eigenvalues are simple. 

Formula (66) follows now from (68). Furthermore, = ^k( n) is an even function for 
odd n > 3, and an odd function for even n> 3 (by Proposition 4.6, see Step 2). Clearly, 
$1 is an even function (e.g. because it is positive in D). Finally, the fact that $2 is an 
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odd function follows by a minor modification of [T^ proof of Theorem 3.4]: by restricting 
attention to either even or odd functions, one obtains that (see Proposition 4.6 for the 
notation) 

the n-th smallest eigenvalue of A'^'^^'^ does not exceed w{{n7i — 7r/2)^/a^); 
the n-th smallest eigenvalue of A°'^'^ does not exceed w{{miY /a^). 



(70) 



In particular, the smallest eigenvalue of A"'^'^ is less then w{'k'^ /o?), so the only candidates 
are Ai and A2 (as shown in the previous step). Since $1 is an even function, we conclude 
that $2 is odd. The lemma is proved. □ 

Remark 4.13. Alternatively, one could prove that $2 is odd as follows. Clearly, a/i2 > 
TT — 7r/8 > 2.748. Hence, if 62 denote s th e right-hand side of (68) for n = 2, then 
62 < 0.4247rii(6)/a, and by Proposition 4.6, there is an eigenvalue Afc(2) of A°'^'^ in the 
interval I2 = (A2 — 62, X2 + ^2)- By using the results of Step 3, one easily proves that 

A2 + £2 < A4 — 64, 

and so k{2) < 4. Since $1 and $3 are even functions, we must have k{2) = 2, and so $2 
is an odd function. 



Remark 4.14. Clearly, a higher-dimensional analogue of (70) holds true: Let D = 
B{0, a) be a ball in R'^ and V{x) = for x E D, V{x) = 00 otherwise. Fix a d-dimensional 
spherical harmonic Y(u) {\u\ = 1) of order / > 0. Suppose that w{^) = + ly^'^ — 1 (or, 
more generally, that w{^) is an operator monotone function), and let Aq = w{—A), where 
A is the Laplace operator in R''. Let (/^^)^ and A^ denote the n-th smallest eigenvalue 
of —A + V{x) and Aq + V{x), respectively, among the eigenvalues corresponding to the 
eigenfunctions of the form /(|a;|)y(a;/|a;|) (therefore, a/x^ is the n-th zero of the Bessel 
function Jd/2+i-i)- Then A^ < w((/i^)^). The proof again is a minor modification of the 
argument used in |l4t proof of Theorem 3.4]. 



Proof of Theorem \l.S\ Recall that by a scaling argument, we only need to prove the 
statement of the theorem for h = c = m = 1; that is, we prove Theorem |L2[. The proof 



relies on Lemma 4.12 By (58), we have afin > 77r/8 > 2.748 for n > 2, and by (57) 



and (55), a/ii > tt/2 — ■&Tr/{2a) > 7r/2 — 1/a > 1.520 when n = 1 and a > 20. By (66), in 

-a/3)e-'^/3 ^^^^ (l + a/3)e- 



both cases. 



-a/3 



|A„ - A„| < 8.610 



< 8.610 



< 8.610 ?ii±^/5)£Z^!< 7.309 

3n7ra 



n7r/2 - 7r/8)a 
(1 + a/3)e~'^/3 



na 



Formulae iBh and (25) for n > 2 or a > 20 are proved. In the remaining case n = 1 and 



< a < 20, we use a more direct approach. By (33) and the definition of A„ (see (41) 
and (|42|), both Ai and Ai are bounded above by w{7r'^ / {Aa'^)) = (7r^/(4a^) + 1)^/^ - 1. 
Hence, 



lAi-Ail < 



4^ 



1 TT^ f IT Tt'^ 

, < mm — , — - 

2a Vvr2 + 4a2 + 2a ~ \2a 8a^ 



Observe that e°-^^/{6 + 2a) is increasing in a > 0. Hence, for a G (0, 6], 



TT 



< 



37Te' 



;i + a/3)e-'^/3 2a 6 + 2a " 6 + 12 



< 3.869. 
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In a similar manner, e'^/'^/(a(3 + a)) is increasing for a > 6, and so for a E [6, 20], 

^ ^ g 323 

(1 + a/3)e-'^/3 8a2 ~ 8a(3 + a) " 8 ■ 20 ■ 23 



It follows that 



|A„ - A„| < 7.309 



;i + a/3)e-"/3 



na 



(71) 



for all n = 1, 2, ... and a > 0. Since En = A„ + 1 and A„ = + 1)^/^ — 1, formula ( [tI] ) 
implies (25). □ 



Proof of Theorem ] Recall that by (55), i)^ = tt /8 + 0(1/ fi) as fi — )■ oo. Furthermore, 
jln > mr/ (2a) — 7r/(8a), so that i?^,^ = 7r/8 + 0(1/??,) as n — t- oo. Hence, 



2a a 



nn 



2a 8a 



O 



0{l/n) here and below corresponds to asymptotic estimates as n ^ oo. It follows that 



mr TT 
2a 8a 



1 + 



By (71), we obtain that 



A. 



mr TT 
2a 8a 



1 + 



as desired (recall that En = A„ + 1). Simplicity of eigenvalues is a part of Lemma 4.12 , □ 



Verification of Remark In this part the general case is studied, without the scaling 
assumption c = h = m = 1. 

(1) We fix h, c, a and n. When m — )■ 0+, the right-hand side of ^ converges to 
8hc/{an). On the other hand, by ([s]) and the estimates (55) of i?^, we have (as in 



the proof of Theorem 1.1 ) /i„ > {n-K/2 — n/8)h/ (mca), i}^^^ = vr/S + 0{m), and 

0(m) 



I <-> Mr: 

mca \ 2 



h /nn TT 
mca \ 2 8 



as m — )■ 0^. It follows that 

^2 /- 



9 /TTT^ 9 ^/ 1\\ / TITT TT ^/ N 

mc V/in + l=mc [fin + 0[^n )) = ~ ~ 8 



Therefore, passing to a limit m — )■ O'^ in (|9]), we conclude that 



En 



he /nn TT 
T vT ~ 8 



< 



he 



an 



(2) Now h, m, a and n are fixed, and we let c — > oo. The right-hand side of (|9| 
converges to exponentially fast. Furthermore, -(9^ = 0{fi) as /i — 0+, and 
by (Is]), /i„ = 0(l/c) as c — )■ oo. Hence, ■(9^^ = 0(l/c), and so 



h /mr 



mca V 2 
as c — 7- oo. This gives 



h 



mca 



nn 



O 



mc^^M+1 = mc" (1 + j2l/2 + 0(/i^)) 



mc^ + 



/nn\- 
2ma2 V vTJ 



O 
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and (10) follows by Higher order (with respect to 1/c) expansions similar 
to (10) can be easily obtained from ( [9| in a similar manner. 
(3) Denote, for brevity, d = fr^mca. By (|8j), ( [54| ) and again ([s]). 



1 frni n /'T-TT 

/in - 3 — - ?9 — 
a \ 2 \2a 



< 



\2aJ 
1 



2 e{~^^) 



na^ 1 + /in 



na 1 + /i^ V 2a 

By -dfx^ < 'dnn/{2a) < min(n/a, vt/S) < {tt /2)n/{d + n). Furthermore, fin > 

n7r/(2a) — 7r/(8a) > 3ri7r/(8a) >n/d. It follows that 



/iri 



1 /nn ^ 
Finally, since w{jj,'^) = {jj,'^ + 1)^/^ 



< 



n 



1 



< 



2n 



a + n + 71^ (a + nY 
1 is convex, by the mean value theorem. 



1 /mi 



a2 V 2 



\2d) 



n-K , 
< 2 — w' 
- 2d 



\2d 



/2a V 



-1/2 



ynn J 



+ 1 



1 /nn 

d vT 



1 /mr 

d vT 



V2a/ 



V2a/ 



Using (2a/(n7r))2 + l > (l/2)(2a/(n7r) + 1)^ > (2/7r2)(a/n + 1)^, we conclude that 



W 



a2 



V2a/ 



< 



TT 



n 



y/2 a + n 



This and (|71|) proves that 
A„ — md^w ' 



1 (n-K 
a \ 2 



2 \2dJ 



/mx\ 
\2d) 



< 



2-n 



n 



V2 {d + ny 



< 



mc 



na 



7.309(1 + a/3)e~''/^ + 4.443 



n + dy 



Formula (11) follows. 



(72) 



□ 



Proof of Proposition 1.5 It suffices to consider the case c = h 



m 



1. The upper 



bound in (12) for n > 1 is just the upper bound in (33). The lower bound for n = 1 is 
trivial. Let, as usual, w(^) = (.^ + 1)^^'^ — 1. Since w{^)/2 > w(^'^/4:), the lower bound 
for n = 2 also follows from the corresponding lower bound of ( p3| ). Finally, let n > 3. As 
in Step 3 of the proof of Lemma 4.12, using (69) we find that 

2^ 



A. 



w 



{n — l)7r' 
2a 



> {w{fLl) -en)-W 

1)ti 



[n 



> 



TT 



2a 
n/a 




0.200nu{b) 



8 Jl + 7iya^ 



Vl + {{n - l)7r)V(2a) 
- 0.200M(a/3) I > 0, 
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and the lower bound in (12) for n > 3 follows. 
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□ 



4.5. Properties of eigenfunctions. We already know that approximations A„ and En 

are indeed close to true eigenvalues A„ and En. Define ^n{x) = 'Pn{x) /\\^Pn\\L^{D)- In 
this subsection we show that also ^n{x) are good approximations to true eigenfunctions 
^n{x). Recall that (^n(x) is defined by (43), with b = a/3, as in the previous subsection. 
In order to prove Propositions [1.6 and L7, we need the following result. 

Proposition 4.15. The sign of ^n can be chosen so that 

< 1,455 ii±^/31^?i^!^v^. 



Proof. Let an,j = (<^n, $i)i2(£,), so that (fn = J2T=i'^n,j^j in L'^{D). Note that a„j = 
when n and j are of different parity. Since we can replace $„ by — if necessary, with 
no loss of generality we assume that an.n > 0. Furthermore, we let Pn = \\'pn\\L^(D)- We 
have 



I C'^njn^n II ^2^£)j llVn 



(73) 



<2\\LPn- an,n^n\\L2 



(D) ■ 



Observe that by ([T2|, (|58| and ([69|, with w(^) = + lY^^ - 1, 

2^ 



An+2 - K>U) 



> 



n + Ijvr 
2a 

nvr/ (2a 



> 



TT 



2a v^l + (n7r)V(2a)2 " Aa"^ ^ I ^ / {2af 
and for n > 3, in a similar manner, 

TT (n-2)7r/(2a) ^ tt^ 



An ~ An-2 > 

Since = an,n-i 



> 



2a + ((n - 2)7r)V(2a)2 ~ 4a2^ 
0, it follows that 



TT 



V(2a) 



IV^n - Oin,n^n\\L2 



(D) 



\j-n\>2 \ 



4a2v^l + 7r2/(2a)2 



/ b-™|>2 



< 



'2aV4a2 + 7r2^ 



L2(D) 



n,J I 



Finally, since /i„a > 3n7r/8 > 1.178n > 1.178, formula (65) gives 



(1 + a/3)e-''/3 (l + a/3)e-'*/3 

< 4.228 ^ ' J,^ < 4.228 ^ ^ ' ' 



Una 



3/2 



1.178n./a 



(74) 
(75) 



Formulae ( 73 ) , ( 74 ) and ( 75 ) combined together yield 



1.178n^/a 



, (l + a/3)V4a2 + 7r2e-'^/3 ^ 
< 1.455 ^ — A/a , 



n 



as desired. 



□ 
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Proof of Proposition L6 The last statement of the proposition was already proved in 



Lemma 4.12 Hence, we only need to prove (15). 



By a scaling argument, it suffices to consider only the case c = h = m = 1. As in the 
proof of Proposition 4.15, we let /3„ = ||'^n||L2(£)). We have 



L2(D) 



'-(D) 



\l\d) ■ 



By Lemma 4.5, as in (67) 

|/3n-«| = 



IV^n|lL2(D) 



<^+4(l + sin^,J/, 

0.383 4 ■ 1.383 ■ 0.217 1.584 
< — — + z < — — . 



(76) 



In a similar manner, by (55) and (61), 
sin '(9 



^" + 4(1 + sin ^^J J^„ < - + < 0.637 + 0.692/i„. 



Hence, 



a - /3„)$„||^,(^) = -a\< ^^^7=^ < min f 0.637 + 0.692/i„ ) . 



By a simple calculation, the minimum is always less than 1.414. Hence 



\L^{D) 



^ 1 . / 1.584 ^ . flM5^ 1.414\ 

< —= mm — : — , 1.414 ) < mm ( — — ] 

a \ Hn 



' / — I ' 
n A/a / 



the last inequality follows from /i„ > 3ri7r/(8a). By the definition of (fn and fn (see (14) 
and (|47])), 



afnix)\ < Ga„, (a + x) + Gu„ (a - x) 



so that by (60) and ( |6l| ) (and 3ra7r/(8a) < fin ^ mT/{2a) in the last inequality), 

1/2 



((0,(»)) 



. , ■0.383- 0.217 ^ ^Q 333 

< 2 mm , — < mm 



< mm 



/in 47r 
0.532v^ 0.887n 



/in 



Finally, by Proposition 4.15, we have 



(l + a/3)V4a2 + 7r2e-'^/3 

|<^n - /9n'^'n||L2(o) < 1-455 ^ — ^/tt . 



n 



By an explicit calculation, the numerator on the right-hand side does not exceed 5.318, 
and a careful estimate gives that it is also less than 35/a. Hence, 



7.738a/S 50.925 



n nwa 
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By collecting all above estimates, we conclude that 

'7.738v^ 50.925\ 
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1.345Va 1.414\ 

+ mm 



n 



a/o J 



+ min 



n ' n^/a J 
0.532Va 0.887r2 



n 



a 



< min 



9.615v^ 1.414 50.925 0.887n 

+ ^ + 



n 



n\/a 



When a < 0.204n, then the first argument of the minimum is smaller than 0.887n/a, 
which is less than the other argument of the minimum. For a > 0.204n (and n > 1), 

1.414 50.925 0.887r2 1.414v^ 50.925v^ 0.887v^ 54.303v^ 



'a Uy/a a 
and the proof is complete. 



v/0.204a 



□ 



Proof of Proposition Again we only consider the case c = h = m = 1. Recall that 
T{t\ X, y) and To(t; x — y) are the kernel functions of the operators T{t) = exp{—tA) and 
To(t) = exp(— tAg). By definition, T(t)<l>„ = e~'^"*$„. By the formula for the Fourier 
transform of TQ{t; x), we have 



r(i;a;,-)|l22p)= / iTit;x,y)ydy< / (Toit; x - y))' dy 



< 



1 

2^ 
1 

TT 



TT 



ds 



TC 



-2t 



2t + 1 



(77) 



27lt 



Let t = 1/(2A„) and /3„ = ||(y9„||i2(£)). We have 



(75 



Since T(t) is a contraction on L^{D), we have (by ( [43| ) and (60)) 

\\T{t)<^n\\L^(D) < \\'fn\\Loo(^D) < sup (g(-x) (o + x) | + g(x) (a - x) |) 
<||FaL..(r)<1 + G^„(0+)< 1.383; 
for the third inequality, note that q{—x) + g(s) = 1 and q{x),q{—x) > 0. Furthermore, 



\\T{t){(3n^n - Vn)\\Loo(^D) = SUP 



T{t;x,y){(3n^niy) -Vn{y))dy 



D 



< sup (||T(t;X, ■)||l2(o) II /3n$n - </'n|L2 



By Proposition |4.15| ([TTj) and the identity {2t + l)/(2t) = 1 + A„ 

1/2 



|T(t)(/3„<l>„-^, 



"/ IIl°o(d) 



< 



2t + 1 
27rt 



, (l + a/3)y4a2 + 7r2e-'^/3 
1.455 ^ ^/a 



n 



< 0.821 



'1 + Xn)a 



n 



1/2 



1 + a/3)V4a2 + 7r2 e''^/^ 
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By (33), 



and therefore 



[1 + Xn)a a + UTT /2 vr 

^ — < — <a + - 

n n 2 



||T(t)(/3„<l>„,-^„)||^^(^) < 0.821 



^/a+V/2{l + a/3)^/4:^ 



n 



Clearly, the numerator in the right-hand side is a bounded function of a > 0. By a careful 
estimate, we find that it is less than 14. We conclude that 



||T(t)(/3„<l>„-<^„)||^^(^) < 



11.494 



n 



Applying the above estimates to ( 78 ) gives 

e-i/2||$ 



1 /11.494 



+ 1.383 



en 



Finally, by ([76]), - a| < 1.584//i„. Together with fin > 37™/(8a) > 1.178n/a, for 
n > 2 this gives 



1.584a / 1.345 
Pn> \ a- > A / 1 Va . 



1.178n 



n 



and so (for n>2) 



l'^n|lL°=(D) — 



'I - 1.345/n 



11.494 



+ 1.383 



n 



(79) 



This already proves that $n(x) are bounded uniformly in x G -D and n = 1,2,.... To 
prove (16), we need to combine (79) with a similar estimate for small n. 



First, we need a different version of (|77|): 

f oo 



\\T{t;x,-)\\l,^^^<- 



'ds + 



1 



TT 



< 



1 



^ 2TXt 



For t = 1/ (2A„), we obtain 
ll^(i^;a;,-)|li2(B) 



< 




TT TT 




TT 



1 + VAJ < 



2A„(1 + A„) 



TT 



Let m;(0 = (^ + 1)^/^ - 1. Since w{i^){\ + w{i^)) = ^ 1 - /TTF < ^^ the upper 
bound of (33) gives 



2^J 



1/2 



||T(«;x,.)lll.,„<^-=.'^^j j(l + ' 
Arguing as in the previous case, we obtain 

e"^^^ ll'^'nllioofz)) = ll^W'^nllioo/^) = sup / T{t;x,y)^n{y)dy 

xeD Jd 



2 mi ra n 



vr 2a V 2 a 



< sup ( ||T(t;a;,-)||L2(B) ||$n||L2(z)) 

a;6-D ^ 



< 



TT 77, 

2 a 



1/2 
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Hence, 



J- /TT 

l'^n|lL°°{z)) — ve 1^— 



1/4 \n in 
< 1.846 W- 
a y a 



A combination of Q (for n > 16) and Q (for n < 15) gives ||$n||L-(D) < T.SSSa^^/^ 
for all > 1. □ 
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